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The  kinetics  of  the  important  stratospheric  reaction  0(^P)  +  HO2  -* 
OH  +  O2  (1)  has  been  studied  at  298  K  as  a  function  of  N2  diluent  gas 
pressure  from  10  to  500  Torr  and  also  as  a  function  of  temperature  from 
266  to  391  K  in  80  Torr  of  N2.  Pulsed  248.5  nm  KrF  laser  photolysis  of 
O3/H2O2/N2  mixtures  produced  OC^D)  and  OH.  The  Oi^D)  was  rapidly 
quenched  to  0(^P)  by  N2  and  HO2  was  produced  by  reaction  of  OH  with 
H202.  0(^P)  was  monitored  by  time-resolved  resonance  fluorescence.  The 
concentration  of  the  excess  species,  HO2,  was  calculated  from  experiment- 
ally measured  parameters.  This  is  the  first  study  of  this  reaction  under 
pseudo-first  order  conditions  at  pressures  greater  than  a  few  Torr  and  the 


first  time  the  described  techniques  have  been  employed  for  studying 
radical-radical  reactions. 

The  temperature  and  pressure  dependence  of  another  important 
stratospheric  reaction,  0(^P)  +  CIO  ->  CI  +  O2  (2),  has  been  studied  over 
the  ranges  231  to  367  K  and  25  to  500  Torr.  Pulsed  351  nm  XeF  laser 
photolysis  of  CI2/O3/N2  mixtures  produced  CI  in  excess  over  O3.  After  a 
delay  sufficient  for  the  reaction  CI  +  O3  -»  CIO  +  O2  to  go  to  completion, 
a  small  fraction  of  the  CIO  was  photolyzed  by  a  266  nm  Nd:YAG  laser 
pulse.  The  decay  of  0(^P)  in  an  excess  of  a  known  concentration  of  CIO 
was  followed  by  resonance  fluorescence  spectroscopy.  A  few  measure- 
ments of  the  rate  coefficient  for  the  reaction  CIO  +  CIO  -*■  products  were 
also  performed. 

The  rate  coefficients  for  both  reactions  (1)  and  (2)  were  found  to  be 
independent  of  pressure.  Our  results  are  described  by  the  following 
Arrhenius  expressions  in  units  of  10"^^  cm^  molecule'^  s'^:  k^CT)  =  (2.91 
+  0.70)exp[(228  ±  75)/T]  and  k2(T)  =  (1.55  +  0.33)exp[(263  +  60)/T]. 
The  errors  are  2a  and  represent  precision  only.  The  results  are  compared 
with  previous  measurements  and  the  implications  of  these  findings  to 
atmospheric  chemistry  are  discussed. 


VI 


CHAPTER  I 
INTRODUCTION 

The  kinetics  of  gas  phase  reactions  between  free  radicals,  species 
containing  one  or  more  unpaired  electrons,  are  of  importance  from  both 
basic  and  practical  standpoints.  Radical-radical  reactions  often  proceed 
on  an  attractive  potential  energy  surface  with  no  barrier  to  the  formation 
of  a  bound  intermediate.  These  reactions  therefore  offer  relatively  simple 
systems  that  have,  over  the  past  two  decades,  become  amenable  to 
theoretical  studies.  An  in-depth  discussion  of  the  theoretical  basis  for  the 
observed  kinetics  of  radical-radical  reactions  is  beyond  the  scope  of  this 
work  but  can  be  found  elsewhere.^  From  a  more  applied  point  of  view, 
radical-radical  reactions  play  critical  roles  in  complex  chemical  systems 
such  as  in  combustion  and  in  the  chemistry  of  the  Earth's  atmosphere. 

Experimentally,  the  kinetics  of  radical-radical  reactions  are  difficult 
to  study.  Radical-radical  reactions  are  usually  fast  and  in  order  to  study 
the  kinetics  it  is  necessary  to  generate  and  quantify  two  reactive 
intermediates.  Note  that  under  the  definition  given  here  for  free  radicals 
there  are  a  few  species,  such  as  NO2,  NO  and  O2,  that  do  not  have  to  be 
generated  in  situ.     It  is  often  the  case  that  the  self-reactions  of  the 
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radicals,  as  well  as  the  reactions  of  the  radicals  with  reaction  products, 

are  rapid,  making  it  more  difficult  to  quantify  radical  concentrations.  The 
most  commonly  employed  techniques  for  directly  studying  the  kinetics  of 
radical-radical  reactions  are  variations  of  flow-tube  experiments.  In  flow- 
tube  studies,  the  concentration  of  one  species  is  monitored  at  a  known 
distance  downstream  from  where  it  is  mixed  with  an  excess  of  a  second 
species.  Knowledge  of  the  flow  velocity  allows  the  distance  to  be 
converted  to  a  reaction  time  and  variation  of  the  distance  and  the 
concentrations  leads  to  the  kinetic  information.  A  major  advantage  of 
flow-tube  experiments  is  adaptability.  Many  different  detectors  can  be 
used  to  obtain  reactant  concentrations  and  a  wide  variety  of  radicals  can 
be  generated,  often  in  an  electrical  discharge  or  by  thermal  decomposition. 
Standard  discharge-flow  methods  are  restricted  to  pressures  in  the  1  to 
10  Torr  range,  although  some  recent  studies  have  reported  experiments 
performed  at  higher  pressures. ^^  A  major  disadvantage  is  the  possibility 
of  undetected  heterogeneous  effects.  A  more  detailed  description  of  flow- 
tube  experiments  can  be  found  elsewhere."*  In  part,  the  present  studies 
were  undertaken  to  develop  an  alternate  experimental  method  for  the 
study  of  radical-radical  reactions.  Previous  results  on  the  two  reactions 
studied  in  the  present  work  are  derived  mainly  from  flow-tube 
experiments.  Our  adaptation  of  the  laser  flash  photolysis-resonance 
fluorescence  (LFP-RF)  technique  to  the  study  of  radical-radical  reactions 
has  allowed  these  two  reactions  to  be  studied  at  higher  pressures  than 


3 
before  and  without  the  same  systematic  errors  inherent  in  discharge  flow 

studies.  A  discussion  of  the  innovations  in  the  LFP-RF  technique,  as  well 

as  comparison  with  previous  results,  will  appear  in  subsequent  chapters. 

The  two  reactions  studied  in  this  work 

0(3p)  +  HO2  -*  OH  +  O2  (1) 

0(3p)  +  CIO  -  CI  +  O2  (2) 

are  both  of  importance  in  the  earth's  stratosphere.  The  major 
stratospheric  chemistry  problem  relates  to  ozone  which,  along  with 
molecular  oxygen,  forms  the  protective  ultraviolet  absorption  layer  under 
which  life  evolved  on  Earth.  Our  knowledge  of  the  natural  dynamic 
balance  between  solar-initiated  photochemistry,  free  radical  reactions,  and 
the  background  gases  that  are  precursors  of  the  free  radicals  in  the 
stratosphere  can  be  found  described  in  detail  in  many  sources. ^"^^  Much 
of  this  knowledge  has  been  gained  in  the  last  two  decades  from  research 
spurred  on  by  the  recognition  that  the  activities  of  man  are  altering  the 
"natural"  composition  of  the  atmosphere.  Of  particular  importance  is  the 
problem  of  how  the  introduction  of  additional  trace  species  to  the 
stratosphere  changes  the  balance  between  the  formation  and  destruction 
of  ozone. 

The  formation  of  a  roughly  constant  level  of  ozone  in  the  strato- 
sphere was  first  explained  by  Chapman^^  via  the  following  mechanism: 
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02  +  hv  (<  190nm)  -20  (3) 

O  +  O2  +  M  -»  O3  +  M  (4) 

O3  +  hv  (<  310  nm)  -  O  +  O2  (5) 

O  +  O3  -  2  O2  (6) 

where  M  represents  a  third-body  collider  (N2  and  O2).  This  oxygen-only 
mechanism  predicts  concentrations  of  odd  oxygen,  O3  and  O,  that  are  too 
large  when  compared  with  measured  levels.  It  is  now  known  that  other 
major  loss  processes  are  extant  in  the  stratosphere.  These  include  the 
catalytic  cycles 

X  -I-  O3  -  XO  +  O2 
XO  -I-  O  -  X    +  O2 


O  -K  O3  -  2  O2  net 

where  X  =  OH,  H,  NO,  CI  or  Br.  It  should  be  noted  that  the  net  result  of 
this  cycle  is  the  equivalent  of  reaction  (6).  At  lower  altitudes,  where  the 
partitioning  between  "odd"  oxygen  species  (O3  and  O)  is  shifted  toward 
the  molecular  species,  the  efficient  ozone  destruction  cycle  becomes: 
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OH  +  Og  ->  HO2  +  O2  (7) 

HO2  +  O3  -*  OH  +  2  O2  (8) 


2  O3  -  3  O2        net 


The  efficiency  of  each  cycle  toward  odd  oxygen  destruction  is 
dependent  upon  the  rate  coefficient  for  the  rate-determining  step  in  each 
cycle  and  the  coupling  with  reactions  that  remove  the  catalytic  radicals. 
Examples  of  the  latter  include 

CI  +  CH4  -  HCl  +  CH3  (9) 

and  OH  +  HO2  -*  HgO  +  O2  (10) 

The  major  hydrogen-oxygen  free  radicals,  OH  and  HO2,  play  critical 
roles  in  the  chemistry  of  the  stratosphere  either  in  catalytic  cycles  or 
through  reactions  that  partition  other  species  between  free  radicals  and 
chemically  stable  forms.  The  dominant  source  of  HO^  radicals  throughout 
the  stratosphere  is  the  reaction 

H2O  +  0(lD)  ^  2  OH  (11) 
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where  the  excited  oxygen  atoms  are  formed  by  the  photolysis  of  O3. 

Reaction  (1)  is  the  rate-determining  step  in  a  catalytic  cycle  that  leads  to 
O3  destruction 

OH  +  O3  -*  HO2  +  O2  (7) 

O  +  HO2  ^  OH  +  O2  (1) 


O  +  O3  -  2  O2  net 


Figure  1  (taken  from  reference  12)  shows  the  conversion  rate  between  OH 
and  HO2  as  a  function  of  altitude  for  all  major  reactions  thought  to 
interconvert  these  species.  Note  that  reaction  (1)  is  the  dominant 
reaction  converting  HO2  to  OH  above  40  km,  i.e.,  in  the  upper 
stratosphere. 

In  1974,  Stolarski  and  Cicerone^^  suggested  that  chlorine  introduced 
into  the  atmosphere  in  rocket  exhaust  could  destroy  O3  by  another 
catalytic  cycle 

CI  +  O3  ^  CIO  +  O2  (12) 

O  +  CIO  ^  CI  +  O2  (2) 

O  +  O3  -  2  O2  net 
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8 
In  this  cycle,  reaction  (2)  is  rate-determining.  Subsequently,  Molina  and 

Rowland^^  identified  chlorofluorocarbons  (CFCs)  as  another,  much  larger 
anthropogenic  source  of  stratospheric  CI.  In  general,  CFCs  are  defined 
as  hydrocarbons  in  which  all  the  hydrogen  atoms  have  been  replaced  by 
either  chlorine,  fluorine  or  a  combination  of  both.  Since  World  War  II, 
CFCs  have  been  produced  for  use  in  a  variety  of  industrial  applications 
such  as  solvents  and  refrigerants.  Relatively  inert  to  oxidation  or  other 
loss  processes  in  the  troposphere,  the  CFCs  mix  into  the  stratosphere 
where  their  ultraviolet  photolysis  produces  free  chlorine  atoms. 
Approximately  80%  of  the  present  atmospheric  load  of  chlorine  is  from 
CFCs  with  the  remainder  coming  from  naturally  occurring  CH3CI. 

Recent  measurements  have  clearly  shown  that  indeed  the  level  of 
stratospheric  O3  has  been  negatively  perturbed.  Trend  estimates  in  the 
loss  of  ozone  in  the  mid-latitudes  as  measured  by  several  methods  and 
predicted  by  model  calculations  are  shown  in  Figure  2  (taken  from 
reference  15).  The  measurement  techniques  include  ground  based 
(Umkehr),  balloon-borne  (Payerne  and  Hohenpeissenberg),  and  satellite- 
borne  instrumentation  (SAGE).  As  can  be  seen  in  the  figure,  the 
measured  loss  of  O3  is  greatest  in  the  lower  stratosphere  between  15  and 
20  km.  More  dramatic  losses  in  the  ozone  column  have  been  observed  for 
the  past  several  years  over  the  continent  of  Antarctica  in  the  austral 
spring.    Heterogeneous  processes  that  take  place  on  polar  stratospheric 
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10 
clouds  freeing  highly  photoactive  chlorine  species  (HOCl  and  CI2)  followed 

by  reactions  of  the  oxides  of  chlorine  are  thought  to  be  responsible  for  the 

so-called  Antarctic  "ozone  hole." 

The  specific  roles  of  reactions  (1)  and  (2)  in  the  chemistry  of  the 

Earth's  atmosphere  are  discussed  further  in  the  following  chapters.   The 

study  of  the  room  temperature  rate  coefficient  for  the  O  +  HO2  reaction 

is  presented  in  Chapter  II,  followed  by  the  temperature  dependence  study 

of  the  same  reaction  in  Chapter  III.     The  study  of  the  temperature 

dependence  of  the  O  +  CIO  rate  coefficient  is  described  in  Chapter  IV. 

Except  for  minor  editorial  changes  necessary  to  integrate  the  three 

separate  studies  into  a  single  thesis,  the  material  in  each  of  these  chapters 

is    taken,    by    permission,    from    three    separately    published   journal 

articles. ^^"^®     A  concluding  chapter  contains  a  summary  of  the  work 

presented  here  as  well  as  a  discussion  of  pertinent  information  that  has 

appeared  in  the  literature  subsequent  to  the  publication  of  the  three 

studies.  Finally,  there  are  two  appendices.  The  first  contains  additional 

details  of  the  experimental  and  calculational  techniques  and  the  second 

appendix  lists,  in  numerical  order,  the  complete  set  of  reactions. 


CHAPTER  II 

PULSED  LASER  PHOTOLYSIS  STUDY  OF  THE  REACTION 
BETWEEN  0(3p)  AND  HO2 


Introduction 
Bimolecular  reactions  involving  two  free  radicals  are  of  great  interest 
because  both  reactants  have  unpaired  electrons  and  hence  could  interact  at 
distances  longer  than  those  typical  of  radical-molecule  encounters.  Also, 
because  of  the  attractive  nature  of  the  encounter,  the  energy  barriers  for  the 
reaction  can  be  nonexistent  or  small;  i.e.,  the  reaction  can  proceed  on  a  purely 
attractive  potential  surface.  Consequently,  association  reaction  channels  can 
become  possible  even  when  there  exist  other  pathwaj^  such  as  simple  atom 
transfers.  The  majority  of  radical-radical  reactions  have,  until  now,  been 
studied  using  the  low  pressure  discharge  flow  technique'^  which  is  extremely 
well  suited  for  these  types  of  reactions  since  the  two  free  radicals  can  be 
created  in  physically  separate  regions  and  then  mixed  to  observe  the  reaction 
of  interest.  However,  the  pressure  range  of  discharge  flow  studies  is  presently 
limited  to  a  few  Torr.  (Currently,  a  few  laboratories  are  developing  high 
pressure  discharge  flow  apparatus.)  Therefore,  it  is  obvious  that  if  the 
association  pathways  in  radical-radical  reactions  are  to  be  observed,  especially 
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when  they  are  in  competition  with  other  pressure  independent  channels,  these 

reactions  will  have  to  be  studied  at  much  higher  pressures. 

In  addition  to  the  above  mentioned  rationale  for  studying  radical-radical 
reactions  at  high  pressures,  there  is  also  a  more  practical  reason  for  such 
studies.  Most  radical-radical  reactions  that  are  important  in  the  atmosphere 
and  in  combustion  systems  take  place  at  pressures  much  greater  than  a  few 
Torr.  If  one  intends  to  model  such  complex  systems,  it  is  necessary  to  ensure 
that  the  input  rate  constants  are  those  applicable  under  the  pressure  conditions 
actually  encountered.  Also,  one  cannot  ignore  the  possibility  of  surface 
enhanced  reactions  influencing  the  results  of  discharge  flow  studies  even 
though,  in  principal,  experimental  checks  for  such  effects  can  be  made.  Lastly, 
multiple  studies  using  vastly  different  techniques  are  necessary  to  recognize 
and  then  (hopefully)  minimize  systematic  errors.  The  last  point  is  especially 
important  since  a  great  deal  of  reliance  is  placed  on  the  absolute  accuracy  of 
kinetic  data  in  modeling  atmospheric  and  combustion  processes. 

In  the  past,  a  few  high  pressure  studies  of  radical-radical  reactions  have 
been  carried  out.  However,  except  for  the  cases  of  self-reactions,  i.e.,  R  +  R 
-^  products  (where  R-  =  free  radical),  most  of  these  studies  have  not  been  able 
to  isolate  the  reaction  of  interest  and  therefore  had  to  rely  on  modeling  a 
complex  scheme  of  reactions  to  extract  the  rate  coefficient  of  interest.  To 
worsen  the  situation,  in  many  cases  the  modeling  had  to  include  reactions 
which  were  themselves  not  well  understood.  In  some  of  these  studies  the 
modeling  calculations  showed  that  the  experimentally  measured  parameters 
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were  not  even  very  sensitive  to  the  value  of  the  rate  constant  of  interest. 

Therefore,  it  is  clear  that  experimental  techniques  which  can  directly  measure 
radical-radical  reaction  rate  constants  at  high  pressures  are  needed. 

In  our  laboratory,  we  are  developing  a  method  based  on  laser  photolysis 
to  selectively  produce  free  radicals  in  the  homogeneous  gas  phase  in  such  a  way 
as  to  isolate  the  reaction  of  interest  and  subsequently  follow  the  course  of  the 
reaction  using  spectroscopic  techniques.  This  method  has  become  feasible 
because  of  advances  in  laser  technology  as  well  as  improvements  in  oxir 
understanding  of  the  photochemistry  of  precursors  and  kinetics  of  reactions 
from  which  the  radicals  of  interest  are  generated.  The  present  paper  describes 
the  first  of  these  studies  where  the  rate  coefficient  for  the  reaction  of  0(^P) 
with  HO2, 

0(3P)  +  HO2  -*  OH  +  O2  (1) 

has  been  measured  at  N2  pressures  ranging  from  10  to  500  Torr. 

Reaction  (1)  was  chosen  for  the  initial  study  for  the  following  four 
reeisons:  (a)  there  is  only  one  set  of  products  possible  for  this  reaction,  at  least 
at  low  pressures;  (b)  the  reaction  could  be  studied  using  only  one  photolysis 
laser;  (c)  there  have  been  two  recent  studies  of  reaction  (1)  at  low  pressures ^^'^^ 
which  are  in  agreement  (see  below)  and  thus  provide  a  good  basis  for 
comparison;  and  (d)  the  chemistry  in  the  system  could  be  kept  under  very  good 
control.  Based  on  the  above  description,  one  might  get  the  impression  that  this 
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technique  is  of  limited  use.  However,  it  will,  in  time,  become  evident  that  with 

carefully  thought  out  systems  and  further  improvements  in  our  understanding 
of  the  photochemistry  of  important  radical  precursors,  this  technique  will  prove 
to  be  extremely  powerful. 

Reaction  (1)  has  been  studied  at  298  K  by  many  investigators. ^^"^^  The 
measured  values  of  k^  range  from  2.7  to  7.0  x  10'^^  cm^  molecule'^  s"^.  All 
studies,  ^^'^^  except  one,^^  have  been  carried  out  using  discharge  flow  methods 
of  pressures  lower  than  4  Torr.  Some  of  these  investigations^^'^^  were  indirect, 
i.e.,  ki  was  measured  relative  to  other  rate  coefficients.  The  rate  constants  for 
the  reference  reactions  themselves  have  recently  been  revised.^"*  The  high 
pressure  study^^  was  very  indirect  and  involved  the  pulse  radiolysis  of  a 
mixture  of  O2,  H2,  and  Ar  followed  by  extensive  modeling  to  extract  k^.  The 
most  recent  low  pressure  studies  by  Keyser^^  and  Sridharan  et  al.,^^  however, 
were  carried  out  such  that  reaction  (1)  was  isolated  and  kj  measured  under 
pseudo-first  order  conditions  to  be  (6.1  ±  0.4)  x  10'^^  and  (5.4  ±  0.9)  x  10'^^ 
cm^  molecule'^s"^,  respectively.  Keyser  has  also  measured  k^  as  a  function  of 
temperature  and  obtained  a  small  negative  activation  energy  for  reaction  (1). 

In  the  study  described  here,  0(^P)  and  HO2  were  produced  by  cophotolysis 
of  O3  and  H2O2  in  N2  at  248.5  nm  using  a  KrF  excimer  laser.  Reaction  (1)  was 
studied  under  pseudo-first  order  conditions  in  [0(^P)]  with  [HO2]  >  [0(^P)]. 
[HO2]  was  calculated  from  the  knowledge  of  various  parameters.  The  rate 
coefficient  for  reaction  (1)  was  measured  at  seven  different  total  pressures 
between  10-500  Torr  where  the  bulk  of  the  gas  was  N2  and  kj  was  found  to  be 


15 
independent  of  pressure  with  an  average  value  of  (6.2  +.  1.1)  x  10"^^  cm^ 

molecule'^s"^. 

Experimental  Section 

Reaction  (1)  was  investigated  at  high  pressures  (10-500  Torr)  where  the 

reaction  was  isolated  and  was  studied  under  pseudo-first  order  conditions  in 

[0(^P)]  with  [HO2]  in  excess  over  [0(^P)].  The  experimented  approach  used  in 

this  investigation  involved  248.5  nm  KrF  excimer  laser  photolysis  of  a  mixture 

of  H2O2  and  O3  in  N2  (or  Ar)  to  produce  known  concentrations  of  HO2  and 

0(^P).   The  temporal  profile  of  [0(^P)]  was  subsequently  followed  using  time 

resolved  resonance  fluorescence  detection.    The  cophotolysis  of  H2O2  and  O3 

in  excess  N2  at  248.5  nm  produces  HO2  and  0(^P)  via  the  following  reactions: 


hv 

H2O2 >  20H  (13) 

248.5 


O 


3 


hv 

->  0(lD)  +  02(alAg)  (14a) 


248.5 


0(3P)  +  02(X32 -)  (14b) 


OH  +  H2O2  — >  HO2  +  H2O  (15) 


0(lD)  +  N2 >  0(3P)  +  N2  (16) 
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[0(^P)]  decays  due  to  the  following  reactions: 


0(3p)  +  HO2  -  OH  +  O2  (1) 


0(3p)  +  H2O2 


^  OH  +  HO2  (17a) 


— ^  H2O  +  O2  (17b) 

0(3p)  ^  loss  from  detection  zone  due  to  diffusion 

and  reaction  wdth  background  impurities  (18) 


The  concentration  of  HO2  was  not  directly  measured  but  was  calculated 
based  on  experimentally  measured  and  other  known  parameters.  The  criteria 
for  producing  a  known  concentration  of  HO2  are  listed  below  and  it  is  worth 
pointing  out  that  the  accuracy  of  the  measured  value  of  k^  directly  depends  on 
how  well  each  one  of  these  criteria  is  met:  (i)  the  concentration  of  the 
photolytic  precursor  [H2O2]  should  be  measured  in  the  system;  (ii)  the 
absorption  cross  section  of  H2O2  at  the  photolysis  wavelength  a(H202^248.5 
nm)  should  be  known;  (iii)  ^OH'  ^^^  quantum  yield  for  the  production  of  OH, 
the  precursor  of  HO2,  should  be  known  at  the  photolysis  wavelength;  (iv)  F,  the 
laser  fluence  should  be  measured;  and  (v)  the  laser  fluence  should  be  "spatially 
uniform."  If  all  these  criteria  are  met,  since  the  concentration  of  H2O2  was 
such  that  the  system  was  optically  thin  at  248.5  nm,  the  concentration  of  OH 
formed  can  be  calculated  using  the  relation, 

[OHlo  =  F  X  [H2O2]  X  a(H202,  248.5  nm)  x  Oqh  (D 
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The  concentration  of  H2O2  was  directly  measured  using  absorption  of 

202.6  nm  zinc  lamp  radiation  in  a  76  cm  long  cell  which  was  traversed  by  the 

reaction  mixture  after  it  exited  the  reactor.    The  absorption  cross  section  of 

H2O2  at  202.6  nm^'*  is  4.30  x  10"^^  cm^.  The  absorption  cross  section  of  H2O2 

at  248.5  nm^'^  is  8.63  x  10"^^  cm^  and  the  quantum  yield  for  OH  production 

from  H2O2  at  248.5  nm^^^^  is  2.    In  addition,  it  has  been  shown  that  OH  is 

produced  in  the  electronic  and  vibrational  ground  state. ^^   In  the  presence  of 

excess  N2,  the  rotational  and  translational  temperature  of  OH  would  be  relaxed 

to  that  of  the  bath  gas  within  a  microsecond,  a  period  very  short  compared  to 

the  time  scale  of  the  experiments. 

Spatial  uniformity  of  laser  beams  is  not  easy  to  obtain.  The  beam  has  to 

be  reproducibly  uniform  for  each  pulse  and  not  so  merely  on  the  average.   If 

there  is  a  spatial  gradient  (as  in  the  case  of  a  Gaussian  beam)  or  if  there  are 

regions  of  high  fluence  (hot  spots)  then  reaction  (1)  would  proceed  at  different 

rates  in  different  parts  of  the  beam.   The  0(^P)  detection  zone  is  of  finite  size 

(~2  cm^)  and  hence  the  measured  0(^P)  decay  rate  will  be  an  unknown 

weighted  average  of  the  rates  of  reaction  (1)  in  the  various  volumes  in  the 

detection  zone.  Moreover,  since  one  cannot  directly  measure  the  fluence  in  the 

detection  zone  itself,  if  the  beam  were  nonuniform,  then  the  fluence  measured 

outside  the  reactor  could  be  an  erroneous  representation  of  the  fluence  in  the 

detection  zone.  To  overcome  all  of  these  problems  we  used  a  device  known  as 

a  segmented  aperture  optical  integrator  to  make  the  photolysis  beam  spatially 

uniform.    The  operation  of  this  device  and  the  extent  of  uniformity  of  the 
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integrated  beams  has  been  described  in  an  earlier  publication  from  our 

laboratory.2^  The  spatial  uniformity  of  the  beam  used  in  the  present 
investigation  was  better  than  5%  over  its  entire  cross  section.  The  beam  was 
made  to  essentially  fill  the  reactor,  and  the  laser  fluence  was  constant  over  the 
entire  (-20  cm)  length  of  the  reactor. 

The  laser  beam  fluence  was  measured  as  the  beam  exited  the  reactor  using 
an  EG  and  G  photodiode  based  radiometer  capable  of  measuring  individual 
pulses.  Even  though  the  calibration  of  this  instrument  is  traceable  to  NBS 
standards,  we  carried  out  actinometry  experiments  to  check  the  calibration. 
These  actinometry  experiments  are  described  below. 

Actinometry 

Preliminary  experiments  using  conventional  actinometry  techniques  such 
as  the  aqueous  phase  ferrioxalate  method  or  gas  phase  ozone  photolysis 
followed  by  end  product  analysis  showed  very  poor  precision  as  well  as  fluence 
dependence.  Therefore,  it  was  clear  that  such  methods  are  ill  suited  for 
measuring  fluence  of  high  powered,  short  pulse  width  lasers.  To  overcome  this 
problem,  we  employed  a  very  simple  method  where  we  measured  the  time 
resolved  loss  of  O3  due  to  photolysis  under  geometrical  conditions  identical  to 
those  employed  in  the  actual  kinetics  experiments.  The  experimental  set  up  is 
shown  in  Figure  3.  A  5  cm  diameter  cell  equipped  with  2"  quartz  windows  on 
opposite  faces  to  transmit  the  photolysis  beam  and  1"  windows  on  tubes 
attached  to  the  sides  of  the  cell  to  transmit  the  analysis  beam  (253.7  nm)  at  90° 
to  the  photolysis  beam  was  used.     A  mixture  of  O3  in  N2  was  slowly  flowed 
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Figure  3.  A  schematic  diagram  of  the  apparatus  used  for 
actinometry  experiments  involving  ozone  loss 
measurements  in  real  time.  The  gas  mixture  containing 
ozone  and  N2  were  premixed  from  measured  flow  rates 
and  entered  at  the  port  marked  gases  in.  The  pressure 
in  the  system  was  measured  at  the  port  marked  gases 
out.  The  hatched  area  represents  the  248.5  nm 
photolysis  beam  of  square  cross  section  from  the 
segmented  aperture  optical  integrator.  The  dashed  line 
represents  the  253.7  nm  Hg  radiation,  L  is  the  distance 
the  253.7  nm  beam  traverses  in  the  path  of  the 
photolysis  beam  and  L'  is  the  length  over  which  ozone 
absorbs. 
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through  the  cell.    The  2"    faces  were  covered  by  a  precisely  cut  3  cm  x  3  cm 

aperture  to  clearly  define  the  area  of  the  photolysis  laser  beam  which  traversed 
the  cell.  A  beam  of  253.7  nm  radiation  from  a  mercury  pen  ray  lamp  run  by  a 
regulated  dc  power  supply  traversed  the  cell  through  the  1"  windows,  thereby 
crossing  the  3  cm  wide  photolysis  beam.  It  was  passed  through  a 
monochromator  and  detected  by  a  quartz  enveloped  1P28  photomultiplier  tube. 
The  output  of  the  photomultiplier  tube  was  monitored  using  a  signal  averager 
operating  in  the  peak  height  analysis  mode. 

The  transmitted  253.7  nm  light  intensity  Ig  was  measured  with  the  cell 
flushed  with  N2.  Then  a  constant  amount  of  ozone  was  added  and  the 
transmitted  intensity  I'  measured.  The  ratio  I'/Ig  was  used  to  calculate  the 
concentration  of  ozone  flowing  through  the  cell.  The  photolysis  laser  was 
turned  on  and  the  integrated  beam  allowed  to  photolyze  the  ozone  in  the  cell. 
The  area  of  cross  section  of  the  beam  was  adjusted  to  fill  the  3  cm  x  3  cm 
aperture  in  front  of  the  cell,  thereby  photolyzing  O3  in  an  identical  area  inside 
the  cell.  The  signal  averager  was  pretriggered  and  the  transmitted  253.7  nm 
light  intensity  monitored.  The  flow  rate  of  the  O3/N2  mixture  through  the  cell 
was  such  that  the  cell  was  completely  replenished  between  laser  flashes 
(repetition  rate  =0.1  Hz).  The  observed  temporal  behavior  of  the  transmitted 
253.7  nm  light  intensity  is  shown  in  Figure  4.  Depending  on  the  signal  level, 
64  to  500  flashes  were  averaged.  As  seen  in  Figure  4,  synchronous  with  the 
laser  pulse  there  is  a  sudden  jump  in  the  transmitted  intensity  to  I^  followed 
by  a  slower  increase  to  a  maximum  I2.  Following  this  maximum,  the  intensity 
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Figure  4.  The  temporal  profile  of  the  intensity  of  the  253.7  nm 
beam  transmitted  through  the  cell  including  the  photoly- 
sis area.  Iq  is  the  intensity  with  no  ozone  in  the  cell.  I' 
is  the  intensity  after  ozone  has  been  added  and  its 
concentration  has  equilibrated.  Notice  the  discontinui- 
ties in  both  Y  and  X  axes.  The  arrow  marks  the  time  at 
which  the  photolysis  laser  was  fired.  I^  and  I2  are  the 
intensities  immediately  after  photolysis  and  after  all 
chemistry  has  stopped.  The  decay  from  I2  to  I',  in  the 
time  scale  of  0-8  s,  is  due  to  refilling  of  the  cell.  For  the 
next  laser  pulse,  the  sequence  starts  at  time  0.  I^/I' 
yields  the  initial  ozone  loss  and  consequently  the  fluence 
while  I'/Iq  yields  the  concentration  of  ozone  in  the  cell. 
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slowly  decreases  back  to  I'  as  the  contents  of  the  cell  are  swept  out  and  filled 

with  a  fresh  mixture. 

The  interpretation  of  the  above  result  is  straightforward.  The  sudden 
jump  is  due  to  the  photolysis  of  O^^^  to  produce  OC^D)  +  02(alAg)  and  0(3p) 
+  02(X^2  ").  The  primary  O3  loss  is  assumed  to  occur  with  a  quantum  3deld 
of  unity.^^  The  slower  ozone  loss  is  due  to  the  following  chemistry: 


O3  +  hv 


■^  0(lD)  +  02(alA  )  (14a) 


-»  0(3p)  +  OgCX^Sg)  (14b) 


0(lD)  +  N2  -  0(3p)  +  N2  (16) 

0(3?)  +  O3  ^  2  O2  (6) 

02(alAg)  +  03^  0(3p)  +  2O2  (19) 

It  has  been  shown  that  O3  photolysis  at  248.5  nm  leads  to  both  0(^D)  + 
02(alAg)  and  0(3p)  +  02(X32g-).2®  The  ratio  of  the  quantum  yields  of  0(3?) 
and  0(^0),  4)0(3p)/00(^D)  has  been  shown^^  to  be  -0.1.  Therefore,  the  net 
loss  of  O3  due  to  the  secondary  chemistry  relative  to  the  primary  process,  i.e., 
reaction  (14)  should  be  -2.8.  The  value  we  obtained  is  2.9  ±  0.4  which 
confirms  that  the  system  is  behaving  as  expected. 
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The  primary  quantity  of  interest,  of  course,  is  Ii/I'.    This  quantity  is 

directly  related  to  the  photon  fluence  in  the  monitoring  zone  of  the  cell  by  the 
relation 

Indi/I')  =  [Osliost  X  (7(03,253.7  nm)  x  L  (ID 

where  [OsJiost  is  the  decrease  in  [O3]  due  to  photolysis,  ct(03,253.7  nm)  is  the 
O3  absorption  cross  section  at  253.7  nm  taken^^  to  be  1.12  x  lO'^*^  cm^,  and  L 
is  the  width  of  the  photolysis  laser  beam  (3.0  cm).  Measurement  of  [03]io8t 
leads  to  the  photon  fluence  F  via  the  relationship, 


F  = (III) 

a(03,248.5  nm)  x  [O3] 


Simultaneous  with  these  measurements,  the  transmitted  laser  fluence  F' 
(corrected  as  described  below)  was  monitored  by  the  EG  and  G  radiometer 
masked  by  a  known  area  (1.32  cm^).  The  ratio  of  F/F'  was  the  quantity  we 
were  after  since  using  it  the  photon  fluence  in  the  kinetics  experiments  could 
be  directly  calculated. 

Four  runs  with  varjdng  amounts  of  O3  and  fluence  were  carried  out.  The 
average  measured  value  of  F/F'  was  1.04  +.  0.08  where  the  error  is  2a.  (F  was 
corrected  upwards  by  7.6%  for  the  light  reflected  back  into  the  cell  by  the  back 
window  and  F'  corrected  by  8%  for  the  loss  due  to  transmission  through  O3 
and  the  back  window.)  In  all  kinetics  experiments  the  fluence  measured  by  the 
radiometer  was  corrected  to  take  the  actinometry  results  as  the  primary 
method  for  fluence  measurements. 
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The  elegance  of  this  method  is  clearly  shown  by  substituting  relation  II 
into  III  and  replacing  [O3]  by  ln(Io/I')/{CT(03,253.7  nm)  x  L'},  to  yield, 


Indi/Io)  X  L' 

F  = (IV) 

a(03,248.5  nm)  x  L 


II,  Iq,  L',  and  L  are  all  directly  measured  very  accurately.  (It  should  be 
noted  that  the  spread  in  measured  intensity  was  -0.005%).  The  only  quantities 
whose  accuracies  need  to  be  relied  on  are  (7(03,248.5  nm)  and  the  quantum 
jdeld  for  ozone  dissociation.  Both  of  these  quantities  are  undoubtedly  accurate 
to  better  than  5%  each;  this  leads  to  an  accuracy  of  11%  in  measuring  F. 
Combining  this  with  possible  errors  due  to  slight  nonuniformities  in  the  laser 
beam  yields  an  accuracy  of  -12%  for  the  photon  fluence  measurements. 

Kinetic  Measurements 

A  schematic  diagram  of  the  apparatus  used  for  carrying  out  kinetic 
measurements  is  shown  in  Figure  5.  The  reaction  cell  was  20  cm  long,  3.8  cm 
i.d.,  and  made  of  brass.  Its  inside  was  coated  with  FEP  teflon  and  over  coated 
with  halocarbon  wax.  A  mixture  of  H2O2,  O3  and  N2  was  flowed  through  the 
cell.  The  H2O2  content  of  this  mixture  was  measured  via  its  absorption  of 
202.6  nm  Zn^  radiation  from  a  Zn  lamp  in  a  76  cm  long  absorption  cell  through 
which  the  reaction  mixture  flowed  after  exiting  the  reactor.  In  a  previous  study 
in  our  laboratory^^,  it  was  shown  that  H2O2  does  not  decompose  to  a 
measurable  extent  in  this  reactor.  O3  does  not  interfere  with  the  H2O2 
measurement  since  (a)  its  concentration  is  very  low  and  (b)  the  absorption  cross 
section  of  O3  at  202.6  nm  is  quite  small.^^  All  gas  transfer  lines  between  the 
source  mixtures  and  the  reactor  were  made  of  teflon  connected  by  a  few 
stainless  steel  cajon  fittings  to  minimize  H2O2  decomposition.   The  H2O2/N2 
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mixture  was  generated  by  bubbling  N2  through  a  bubbler  containing  H2O2. 
Ninety  percent  pure  H2O2  (from  FMC  corporation)  was  used  after  N2  had  been 
bubbled  through  it  for  24  to  36  h  to  remove  H2O  and  O2. 

The  O3/N2  mixture  was  stored  in  a  12  liter  bulb.  The  concentration  of  O3 
in  this  mixture  was  measured  using  253.7  nm  absorption.  The  O3  was  prepared 
by  passing  99.99%  pure  O2  (Matheson  Gas  Products)  through  an  ozonizer  and 
collected  on  silica  gel  at  195  K.  Before  use  ozone  was  liquified  at  77  K  and  O2 
pumped  out.  N2  (99.9995%)  and  Ar  (99.9995%)  were  both  purchased  from 
Matheson  Gas  Products  and  used  as  supplied.  The  pressure  in  the  system  was 
measured  as  the  mixture  exited  the  reactor  using  a  capacitance  manometer. 
The  O3/N2  mixture  and  pure  N2  were  flowed  through  mass  flow  meters  and 
their  flow  rates  controlled  by  needle  valves.  Since  H2O2  will  decompose  if 
passed  through  flow  meters,  a  measured  amount  of  N2  was  bubbled  through  the 
bubbler  containing  H2O2  and  the  mixture  thus  generated  was  added  to  the 
O3/N2  stream  flowing  into  the  reactor. 

The  reaction  cell  was  equipped  with  four  ports  in  its  middle.  Two  ports 
were  used  to  transmit  the  output  of  an  oxygen  resonance  lamp  which  was 
filtered  by  a  CaF2  window  to  remove  Lyman-a  radiation.  The  resonance 
fluorescence  excited  by  the  resonance  lamp  beam  was  collected  at  90°  by  a 
MgF2  lens  and  focused  onto  the  photocathode  of  an  EMR-542G  solar  blind 
photomultiplier  tube.  The  output  of  the  PM  tube  was  amplified,  discriminated 
against  noise,  and  fed  into  a  signal  averager  operating  in  the  multichannel 
scaling  mode. 

A  typical  kinetics  experiment  consisted  of  flowing  0.3-1  x  10^^  cm'^  of  O3 
in  10-500  Torr  of  N2  through  the  reaction  cell  and  the  absorption  cell.  J,,,  the 
intensity  of  the  202.6  nm  Zn"*"  line  transmitted  through  the  absorption  cell  was 
measured.  Next,  the  H2O2/N2  mixture  was  introduced  into  the  stream  of  O3/N2 
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mixture,  and,  after  equilibration  of  H2O2  with  the  walls,  J,  the  transmitted 
202.6  nm  light  intensity,  was  measured.  Using  the  values  of  J  and  Jq,  [H2O2] 
was  directly  calculated;  it  ranged  from  0.3  to  2  x  10^^  cm'^.  J  was  continuously 
monitored  throughout  the  course  of  the  experiments.  The  signal  averager  was 
pretriggered  10  ms  before  the  photolysis  laser  was  fired  to  obtain  the 
background  count  rate  and  the  delay  was  precisely  controlled  by  a  delay 
generator.  The  photolysis  beam,  which  had  been  made  spatially  uniform  by  the 
optical  integrator,  entered  the  reactor  and  initiated  the  chemistry.  [0(^P)]  was 
monitored  as  a  function  of  time.  Typically,  16  to  512  individual  0(^P)  temporal 
profiles  were  averaged  to  obtain  one  curve.  During  the  averaging  period,  the 
fluence  of  each  pulse  of  the  photolysis  beam  was  measured  after  the  beam 
exited  the  reactor,  using  the  (calibrated)  E  G  and  G  radiometer.  The  laser 
fluence  was  then  varied  to  produce  a  different  concentration  of  HO2.  Typically, 
five  to  eight  fluence  values  were  used  to  create  1-10  x  10^^  cm"^  of  HO2.  After 
this  series  of  runs  were  completed,  the  H2O2  flow  was  turned  off  and  Jq 
remeasured.  It  should  be  noted  that  in  these  series  of  runs,  for  each  value  of 
HO2,  the  ratio  of  [H02]o/[O^P)]o  is  not  altered.  To  change  this  ratio,  the 
composition  of  the  mixture  was  varied. 

During  the  course  of  each  kinetic  run,  it  is  necessary  that  the  fluence  (per 
pulse)  be  constant.  The  KrF  excimer  laser  used  in  the  present  experiments  was 
equipped  with  a  gas  processor  which  continuously  circulates  the  gas  mixture 
through  the  laser  cavity,  cleans  up  the  mixture,  and  replenishes  F2  such  that 
the  laser  energy  is  extremely  stable.  The  measured  peak  to  peak  fluctuations 
were  typically  -2%,  and  the  average  of  10  to  512  pulses  much  more  stable. 
This  stability  made  it  unnecessary  to  include  an  energy  monitoring  device  into 
the  circuitry  triggering  the  signal  averager  to  discard  any  individual  temporal 
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profile  obtained  for  a  given  laser  pulse  which  had  not  met  the  energy 
requirement. 

Results  and  Discussion 
All  experiments  were  carried  out  under  pseudo-first  order  conditions  in 
[0(^P)]  with  [HO2]  in  excess  over  [0(^P)].  As  will  be  pointed  out  later,  in  this 
particular  system  it  is  not  essential  to  keep  [HO2]  much  greater  than  [0(^P)]. 
If  the  postphoto lysis  chemistry  is  governed  only  by  reactions  (15)  -  (18)  and 
reaction  (1),  then  the  temporal  behavior  of  [6(3?)]  (after  all  OH  has  been 
converted  to  HO2)  is  given  in  Equation  (V). 


-d[0(3p)] 

-  =  (kiLHOa]  +  k'd)[0(3p)]  (V) 


dt 


where  k'd  =  ki7[H202]  +  k^s-  For  constant  concentrations  of  H2O2  and  HO2, 
Equation  (V)  can  be  integrated  to  yield 

[0(3P)]  =  [0(3p)]oexp{(-ki[H02]  +  k'd)t}  (VI) 

Figure  6  shows  plots  of  ln[0(^P)]  versus  time  for  three  different  values  of 
[HO2]  at  fixed  values  of  [H2O2]  and  [O3].  It  is  clear  that  Equation  (VI)  is 
obeyed.  The  slope  of  these  lines  3deld  k'  =  ki[H02]  +  k'd-  Figure  7  shows  a 
plot  of  k'  versus  [HO2]  and  the  slope  of  this  line  yields  kj.  The  value  of  kj 
measured  using  N2  as  the  diluent  gas  under  various  experimental  conditions 
are  listed  in  Table  I.  Each  of  these  values  were  derived  by  measuring  at  least 
five  different  k'  values. 
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Figure  6.  Temporal  profile  of  0(^P)  following  the  photolysis  of  a 
mixture  of  H2O2,  O3  and  100  Torr  N2  at  298  K.  The 
concentrations  of  HO2  were  varied  by  changing  the 
fluence,  and  are  shown  next  to  each  curve.  For  the  sake 
of  visual  clarity,  the  curves  are  displaced  on  the  Y  axis. 
Different  numbers  of  traces  were  averaged  for  each 
curve.  The  slope  of  each  line  gives  k',  the  measured 
pseudo-first  order  rate  constant. 
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CHO2],10^2  cm 


-3 


Figure  7.  Plot  of  k'  (the  measured  pseudo-first  order  rate 
constant)  versus  [HO2]  at  100  Torr  N2  pressure  and  298 
K.  Each  [HO2]  has  been  corrected  for  the  occurrence  of 
the  OH  +  HO2  reaction  as  described  in  the  text.  The 
slope  of  the  line  yields  k^.  The  intercept  represents  k^ 
=  ki8  +  ki7  [H2O2]. 
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One  series  of  experiments  was  carried  out  with  100  Torn  of  Ar  as  the 
diluent  gas.  In  these  experiments,  a  large  fraction  of  the  0(^D)  produced  by  O3 
photolysis  reacted  with  H2O2  to  produce  HO2, 

0(lD)  +  H2O2  -  OH  +  HO2  (20) 

OH  +  H2O2  -  H2O  +  HO2  (15) 

This  HO2  production  could  not  be  helped  since  argon  is  a  very  poor 
quencher^^'^^  of  0(^D)  and  hence  reaction  (20)  was  competitive  with  the 
quenching.  The  extra  production  of  HO2  (2  H02's  for  each  0(^D))  accounted 
for  -12%  of  the  net  HO2  produced.  The  value  of  kj  obtained  in  100  Torr  of  Ar 
(after  correcting  for  certain  reactions  affecting  [H02]o  which  will  be  discussed 
later)  is  (6.27  +  0.68)  x  lO'll  cm^  molecule"!  s"!. 

It  is  necessary  to  mention  that  some  preliminary  measurements  of  k^  were 
carried  out  using  266  nm  photolysis  (4th  harmonic  Nd:Yag)  in  75  Torr  of  N2. 
These  experiments  yielded  a  value  of  (7.96  +.  0.65)  x  10"^^  cm^  molecule'^s'^, 
(errors  are  2a  and  do  not  include  systematic  errors)  which  were  reported  at  a 
recent  meeting.^  Since  then  we  have  discovered  that  two  corrections  have  to 
be  made  to  this  number,  both  due  to  the  measured  fluence.  First,  the 
actinometry  experiments  decrease  this  value  of  kj  by  4%.  Second,  we  have 
neglected  to  correct  for  the  loss  in  measured  fluence  due  to  the  back  window 
as  well  as  the  increase  of  the  fluence  in  the  reaction  zone  due  to  the  reflection 
off  the  back  window;  this  further  decreases  the  value  by  -16%,  yielding  k^  = 
(6.5  ±  0.6)  X  10'!!  cm^  molecule'^s'^  in  reasonably  good  agreement  with  our 
other  results. 
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Secondary  Reactions  that  Effect  FHOgl 

As  mentioned  earlier,  during  the  course  of  this  work  [HO2]  was  not 
directly  measured  but  calculated.  Therefore,  a  great  deal  of  attention  was  paid 
toward  recognizing,  eliminating,  and  when  unable  to  eliminate,  accounting  for 
HO2  production  or  depletion  reactions.  Each  of  the  reactions  that  can  alter  the 
calculated  [HO2]  is  considered  and  its  possible  effects  on  the  measured  value 
of  ki  are  discussed  in  this  section. 

One  of  the  nice  features  of  the  chemical  system  used  in  this  study  is  that 
HO2  should  be  regenerated  after  consumption  in  reaction  (1): 

0(3P)  +  HO2  -  OH  +  O2  (1) 

OH  +  H2O2  -  HO2  +  H2O  (15) 

To  check  the  effectiveness  of  the  regeneration,  a  few  experiments  were  carried 
out  with  [H02]o/[0(^P)]o  ratio  as  low  as  3.  Indeed,  the  0(^P)  decay  plots  were 
exponential  and  the  value  of  k^  thus  obtained  was,  within  experimental  error, 
identical  to  those  measured  at  higher  values  of  this  ratio  (see  Table  I). 
However,  most  of  the  experiments  were  carried  out  with  [H02]o/[0(^P)]o 
greater  than  10,  with  an  average  value  of  15. 

Even  though  at  first  glance,  looking  at  the  exponential  decays  of  0(^P) 
shown  in  Figure  6  and  the  linear  variation  of  k'  with  [H02]o  shown  in  Figure 
7,  it  seems  that  the  system  is  well  behaved,  a  detailed  consideration  of  the 
chemistry  in  the  system  shows  that  there  are  small  contributions  by  some 
secondary  reactions  that  do  effect  the  calculated  HO2  concentration  and  they 
cannot  be  ignored.     The  contributions  of  these  secondary  reactions  were 
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calculated  and  the  measured  value  of  k^  corrected  for  these  complications.  The 
corrections  due  to  the  secondary  reactions  are  individually  discussed  below. 

OH  +  HO2  reaction.  During  the  course  of  the  reaction,  OH  was  generated 
in  reaction  (1)  and  consumed  by  reaction  (15).  The  steady  state  concentration 
of  OH  in  the  system  was  typically  2%  of  the  [0(^P)]  and  therefore  no  error  was 
introduced  due  to  the  occurrence  of  the  OH  +  HO2  reaction, 


M 
OH  +  HO2  -  H2O  +  O2  (21) 


However,  immediately  after  the  creation  of  OH  by  H2O2  photolysis  the 
concentration  of  OH  free  radicals  was  very  high  and  the  contribution  from 
reaction  (21)  cannot  be  ignored.  The  consequence  of  reaction  (21)  is  that  some 
of  the  OH  reacts  with  HO2  and  not  with  H2O2,  thereby  leading  to  an  [H02]o 
value  which  is  less  than  that  calculated  by  neglecting  reaction  (21).  To  take 
this  error  into  account,  a  numerical  integration  of  the  kinetics  in  the  post- 
photolysis  chemistry  was  carried  out  and  the  [HO2]  that  is  actually  produced 
was  thus  calculated.  The  integration  technique  is  discussed  in  the  final  section 
of  Appendix  A.  The  following  two  reaction  scheme  was  used  in  the  calculation: 

OH  +  H2O2  -"  HO2  +  H2O  (15) 


M 
OH  +  HO2  -*  H2O  +  O2  (21) 


The  recommended  value^'*  of  k2i, 

k2i  =  (7  +  4  Patni)  X  10-11  cm^  molecule-^s-l 
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where  Patm  i^  pressure  in  atmospheres,  was  used.  The  integrations  showed 
that  [H02]o  calculated  assuming  reaction  (21)  to  be  negligible  was  as  much  as 
8%  higher  than  that  actually  produced  in  some  of  our  experiments.  To  check 
for  this  effect,  one  series  of  experiments  were  carried  out  where  -5  x  10^^  cm"^ 
of  H2O2  was  photolyzed  at  various  fluence  values  to  produce  different  ratios  of 
[H202]/[OH]o,  some  as  low  as  800.  Figure  8  shows  a  plot  of  measured  k'  versus 
cedculated  [H02]o  (open  circles).  It  is  clear  that  the  plot  deviates  from  a 
straight  line  at  higher  [HO2]  (i.e.,  lower  [H202]/[OH]o).  However,  when 
allowance  is  made  for  reaction  (21),  as  estimated  by  the  integration  scheme,  the 
plot  is  linear.  Each  value  of  [HO2]  (in  the  kinetic  runs)  was  corrected  to  allow 
for  the  occurrence  of  reaction  (21).  The  corrections  to  [HO2]  varied  from  0.5% 
to  8%  and  was  on  the  average  -3%.  It  should  be  noted  that  the  corrections 
were  larger  for  higher  [HO2]  since  it  was  at  these  values  that  the  [H202]/[OH]o 
ratios  were  lowest.  The  contribution  due  to  reaction  (21)  could  be  made 
negligible  by  making  the  H2O2  concentration  very  large.  However,  three 
factors  discouraged  such  increases.  First,  H2O2  absorbs  130.3  nm  oxygen 
resonance  radiation  and  increasing  [H2O2]  beyond  ~2  x  10^^  cm'^  resulted  in 
prohibitively  degraded  sensitivity  for  oxygen  atom  detection.  Second,  the  loss 
of  0(^P)  due  to  reaction  with  H2O2  (reaction  17)  would  also  increase,  thereby 
making  k'^  large.  Third,  the  fraction  of  0(^D)  reacting  with  H2O2  relative  to 
that  being  quenched  by  N2  would  increase.  Therefore,  [H2O2]  was  always  held 
at  less  than  2  x  10^^  cm'^. 

0(^D)  +  H2O2  reaction.  O3  photolysis  at  248.5  nm  produces  0(^0)  90% 
of  the  time.^^  In  the  presence  of  a  large  concentration  of  N2  the  majority  of 
0(^D)  is  quenched  to  0(^P).  However,  a  small  fraction  of  0(^D)  can  react  with 
H2O2  to  produce  two  molecules  of  HO2  for  each  0(^D)  undergoing  such  a 
reaction.      The  extent  to  which  reaction  (20)  can  contribute  towards  the  net 
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Figure  8.  Plot  of  k'  versus  [HO2]  at  250  Torr  and  298  K.  The 
open  circles  represent  [HO2]  uncorrected  for  the 
occurrence  of  the  OH  +  HO2  reaction.  The  filled  circles 
were  obtained  by  correcting  the  HO2  concentration  as 
described  in  the  text.   The  slope  of  the  line  yields  kj. 
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HO2  produced  depends  on  the  ratio  [H202]/[N2]  as  well  as  the  ratio 
[03]/[H202].  In  all  experiments  where  N2  was  the  diluent  gas,  except  those 
with  10  Torr  of  N2,  [HO2]  produced  via  the  reaction  sequence  (20)  followed  by 
reaction  (15)  could  be  made  less  than  0.05%  of  that  produced  via  H2O2 
photolysis.  In  the  10  Torr  experiments,  however,  the  [HO2]  produced  by  0(^D) 
reaction  ranged  from  0.6%  to  2%  of  the  total  and  this  contribution  was  taken 
into  account.  The  worst  situation  was  in  the  case  of  one  series  of  experiments 
where  argon  was  used  as  the  diluent  gas  as  discussed  earlier.  This  reaction  can 
be  made  negligible  by  working  with  high  [H02]/[0(^P)]o  as  was  the  case  in 
most  experiments. 

0(3p)  -t-  H2O2  reaction.   The  reaction  of  0(3p)  with  H2O2, 

0(3P)  +  H2O2         >  OH  +  HO2  (17a) 

>  H2O  +  O2  (17b) 

can  also  increase  the  [HO2]  during  the  course  of  reaction  (1).  However,  this 
reaction  is  very  slow,  k^  =  -1.45  x  10'^^  cm^  molecule"^s'^  ^^;  so  within  the 
time  0(^P)  decays  to  10%  of  its  initial  value,  the  maximum  amount  of  HO2 
produced  would  be  less  than  -2%  in  the  worst  case,  i.e.,  higher  [H2O2]  (~2  x 
10^^  cm"^)  and  the  lowest  measured  decay  rates  (-70  s'^).  Therefore,  the 
contribution  due  to  this  reaction  was  neglected. 

HO2  reactions  with  impurities.  HO2  is  a  reactive  free  radical  and  hence 
can  react  with  impurities  in  the  system  thereby  leading  to  an  unrecognized 
depletion  of  its  concentration.  Keeping  in  mind  that  the  only  diagnostic  at  our 
disposal  is  the  measurement  of  ©("^P)  temporal  profiles,  in  order  to  check  for 
loss  of  HO2  because  of  reactions  with  impurities  in  the  system,  we  varied  the 
[H2O2]  while  creating  the  same  initial  concentration  of  HO2  by  correspondingly 


38 

varying  the  fluence.  The  measured  0(^P)  decay  rates  were  not  affected  thereby 
indicating  the  lack  of  impurities  coming  in  with  H2O2.  The  HO2  +  O3  reaction 
is  too  slow  to  deplete  [HO2]  and  hence  can  be  ignored.  Variations  in  [O3]  had 
no  effect  on  measured  0(^P)  decay  rates  which  again  indicated  the 
unimportance  of  HO2  reaction  with  O3  as  well  as  with  any  impurities  coming 
in  with  O3.  To  check  for  impurities  being  swept  in  with  the  diluent  gas  N2  the 
flow  rate  of  the  gas  mixture  through  the  cell  was  varied,  and  found  to  have  no 
effect.  This  test,  of  course,  does  not  eliminate  the  presence  of  impurities  in  N2 
itself  We  used  N2  of  the  highest  available  purity  and  it  seems  unlikely  to 
contain  any  reactive  impurities.  In  addition,  the  measured  value  of  kj  was 
independent  of  N2  pressure  from  10  to  500  Torr.  (Of  course,  it  is  possible  that 
there  can  be  an  accidental  cancellation  if  k^  increases  with  pressure  and 
compensates  for  decreases  in  [HO2]  due  to  reaction  of  impurities.  This 
possibility  seems  highly  unlikely.) 

HO2  +  HO2  reaction.  During  the  course  of  the  kinetic  measurements, 
many  observations  indicated  the  importance  of  the  HO2  +  HO2  reaction.  When 
the  0(^P)  detection  sensitivity  was  high  (i.e.,  when  all  windows  were  clean  and 
the  resonance  lamp  tuned  for  minimal  line  reversal),  it  was  clear  that  0(^P) 
decays  were  getting  slightly  nonexponential  at  longer  reaction  times  with  the 
decays  getting  slower  as  reaction  (1)  proceeded.  The  observed  reaction  time 
where  this  departure  was  evident  always  occurred  after  two  1/e  times  and  was 
independent  of  [H02]o-  Also,  these  deviations  were  slightly  more  pronounced 
at  higher  pressures  of  N2  than  at  lower  pressures.  Such  a  decrease  in  decay 
rates  of  [0(^P)]  can  be  due  to  either  the  HO2  concentration  decreasing  as 
reaction  (1)  proceeds  or  due  to  generation  of  0(^P)  via  secondary  reactions.  As 
will  be  discussed  later,  the  second  possibility  can  be  ignored. 
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To  take  HO2  loss  via  the  self-reaction  into  account  we  derive  the  following 
expression  for  the  0(^P)  temporal  profile  based  on  reactions  (1),  (17),  (18),  and 
(22), 

0(3p)  +  HO2  ^  HO  +  O2  (1) 

HO2  +  HO2  -*  H2O2  +  O2  (22) 


0(^P) >  loss  due  to  diffusion  and  reaction  (17). 


[0(3?)]       ki 

In  — = ln(l  +  2k22[H02]ot)  +  kd't  (VII) 

[0(3p)]o      2k22 


The  logarithm  can  be  expanded  into  a  series  since  -  1  <  2k22[H02]ot  <  1,  and 
the  cubic  and  higher  terms  neglected  to  obtain: 


[0(3P)] 
In  -- —  =  kiLHOalot  x  (1  -  k22[H02]ot)  +  k'dt  (VIII) 

[0(3p)]o 


If  [HO2]  is  indeed  time  independent  and  equal  to  [H02]o,  then  a  plot  of 
ln[0(3p)]  versus  time  should  be  a  straight  line  and  its  slope  should  equal 
(kiCHOglo  +  kd').  In  Figure  9,  we  have  plotted  ln([0(3p)]/[0(3p)]o)  versus  time 
(curve  a,  open  circles)  for  one  experiment  at  250  Torr  N2.  It  is  seen  that  up  to 
two  1/e  times  the  decays  do  look  exponential  as  does  the  data  shown  in  Figure 
6.  However,  at  longer  reaction  times,  the  curve  deviates  from  a  straight  line. 
When  each  [0(3p)]j  is  corrected  for  the  depletion  of  [HO2]  using  expression 
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Figure  9. 


(a)  Plot  of  ln[0]/[0]o  versus  reaction  time  (open  circles). 
The  scale  is  on  the  right  as  shown  by  the  arrow.  The 
solid  line  is  the  least  squares  fit  of  this  data  up  to  two 
1/e  times  (i.e.,  down  to  ln[0]/[0]o  =  -2).  The  slope  of 
this  line  is  (424  ±_  20)s"^  where  the  error  is  2a.  Notice 
the  deviation  of  the  points  from  this  line  at  reaction 
times  longer  than  ~6  ms,  (b)  Plot  of  (ln[0]/[0]o)/(l- 
k22[H02]t)  versus  reaction  time  t.  This  correction 
accounts  for  the  loss  in  HO2  due  to  its  self-reaction 
during  the  course  of  reaction  (1).  The  scale  is  on  the 
left.  The  solid  line  is  the  least  squares  fit  of  the  data  up 
to  10  ms.  The  slope  of  the  line  is  (453  ±_  13)  s'^  where 
the  error  is  2a.  Notice  the  linearity  of  the  points.  The 
curve  a  has  been  shifted  up  by  1  ms  for  the  sake  of 
visual  clarity. 
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VIII  and  we  plot  ln([0(3p)]/[0(3p)]o)/(l  -  k22[H02]ot)  versus  time,  the 
calculated  points  indeed  fall  on  a  straight  line  (curve  b).  In  addition,  it  is  seen 
that  the  slope  of  line  b  is  6.8%  larger  than  that  of  curve  a  calculated  using 
points  up  to  two  1/e  times.  It  is  worth  noting  that  curve  b  shows  what  the 
[0(^P)]  decay  would  be  if  HO2  was  not  continuously  decreasing  due  to  reaction 
(22),  but  stayed  constant  at  [H02]o.  The  increase  in  the  calculated  value  of  k'l 
using  the  corrected  decay  curves  over  the  uncorrected  curve  is  -7.8%,  i.e., 
(^'uncorrected  '  k'd)/(k' corrected  "  k'd)  =  1078.  This  is  exactly  what  we  would  get 
if  we  used  the  [HO2]  at  the  point  where  reaction  (1)  has  proceeded  to  its  1/e 
point.  Since  in  all  kinetic  measurements  we  used  the  0(^P)  decays  up  to  two 
1/e  times  in  calculating  k'j  values,  we  can  simply  correct  the  [H02]o  value  to 
represent  that  at  the  1/e  point.  The  value  [HO2]  at  the  1/e  point  for  reaction 
(1),  [H02]'o,  is  related  to  [H02]o  by  the  simple  relation 

[H02]'o  =  [HO]o/{l  -  (2k22A:i)}  (IX) 

Most  of  our  kinetic  data  were  corrected  using  this  method  rather  than 
correcting  each  decay  curve  point  by  point.  One  set  of  data  was  analyzed  by 
correcting  each  decay  curve  and  the  obtained  value  of  k^  was  indistinguishable 
from  that  obtained  using  the  above  simpler  method.  The  correction  of  k^  for 
this  HO2  loss  process  ranged  from  5.4%  at  10  Torr  to  9.7%  at  500  Torr.  Both 
the  uncorrected  (column  5)  and  corrected  values  (column  6)  of  kj  are  listed  in 
Table  I. 

Loss  of  HO2  due  to  diffusion.   It  is  possible  for  HO2  to  diffuse  from  the 
detection  zone  to  the  walls  of  the  reactor  and  be  lost  during  the  course  of  the 
reaction.  To  minimize  this  possibility,  the  photolysis  beam  was  made  to  nearly 
fill  the  reactor.   [0(^P)]  was  monitored  at  the  center  of  the  beam  so  that  HO2 
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loss  due  to  diffusion  would  be  minimal.  Therefore,  HO2  loss  due  to  diffusion 
had  to  be  quite  rapid  to  be  important.  It  seems  unlikely  to  be  faster  than  that 
of  0(^P)  which  was  observed  to  be  less  than  ~5  s'^  in  this  system.  In  addition, 
if  this  loss  was  important,  at  low  [HO2]  the  observed  0(^P)  temporal  profiles 
would  be  nonexponential,  and  be  unobservable  at  high  [HO2].  As  discvissed 
earlier,  slightly  nonexponential  decays  were  observed  at  all  [HO2]  and  were 
completely  accounted  for  by  reaction  (22).  Therefore,  we  conclude  that  HO2 
loss  due  to  diffusion  was  unimportant  and  [HO2]  was  not  corrected  to  take  this 
process  into  consideration. 

Secondarv  reactions  involving  0(^P).  If  0(^P)  is  either  removed  by 
reactions  other  than  (1)  and  (17)  or  created  during  the  course  of  reaction  (1), 
our  measured  \alue  of  k^  would  be  erroneous.  The  reaction  of  0(^P)  with 
H2O2  [reaction  (17)]  does  occur  and,  as  pointed  out  earlier,  merely  contributes 
to  an  increased  intercept.  The  concentration  of  O3  used  in  these  experiments 
(typically  <  1  x  lO^^  cm'^)  cannot  deplete  0(3p)  via  the  0(3p)  +  O3  -  2O2 
reaction.  Just  after  photolysis,  a  small  fraction  of  0(^P)  could  react  with  OH, 

0(3p)  +  OH  -*  O2  +  H  (23) 

However,  the  amount  of  OH  thus  depleted  is  negligible  (-0.05%)  and 
inconsequential  since  the  0(^P)  temporal  profiles  for  the  first  100  ^s  after  the 
laser  pulse  was  not  followed.  Also,  any  0(^P)  lost  by  reaction  (23)  will 
regenerate  an  equivalent  amount  of  HO2  due  to  subsequent  reactions, 

H  +  HO2  -►  2  OH  (H2  +  02,0  +  H2O)20  (24) 

OH  +  H2O2  -*  HO2  +  H2O  (15) 
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Reaction  (23)  cannot  be  an  important  contributor  to  the  0(^P)  loss  during 
reaction  (1)  since  the  steady  state  concentration  of  OH  is  very  low  (-10^  -  10^ 
cm'^)  as  discussed  earlier.  Also,  since  k23  =  ^i,'^^  reaction  of  0(^P)  with  OH 
would  be  indistinguishable  from  that  with  HO2.  The  only  reaction  that  can 
produce  0(^P)  during  our  experiments  is  that  due  to  02(a^Ag)  +  O3  -*■  0(^P) 
+  2O2.  02(a^Ag)  is  formed  by  ozone  photolysis.  The  rate  coefficient  for  this 
reaction^'*  is  too  slow  to  be  important  under  our  experimental  conditions  which 
employed  only  1  x  10^^  cm"^  of  ozone.  Also,  variations  of  [O3]  had  no  effect  on 
the  measured  value  of  kj  implying  the  negligible  importance  of  02(a^Ag). 

Possible  effects  of  H2O.  Reactions  of  HO2  are  known  to  be  enhanced  by 
the  presence  of  H20.^'^  In  our  experiments,  H2O2  was  the  source  of  HO2  and 
it  is  hard  to  remove  water  from  the  system  since  H2O2  decomposition  leads  to 
H2O  formation.  Therefore,  to  check  for  possible  H2O  effects,  two  series  of 
experiments  using  100  Torr  N2  were  carried  out  where  up  to  1.5  Torr  of  H2O 
was  added.  The  measured  values  of  kj  were,  with  the  combined  experimental 
error,  unaltered  and  did  not  show  any  trends  with  H2O  concentration.  We 
could  not  use  higher  concentrations  of  H2O  since  it  absorbs  130.3  nm  oxygen 
resonance  radiation,  thereby  severely  decreasing  0(^P)  detection  sensitivity. 
These  experiments  did  show,  however,  that  our  measured  values  of  kj  could  not 
have  been  affected  by  the  presence  of  H2O  since  the  maximum  amount  of  H2O 
that  could  have  been  in  the  system  was  -10%  of  the  H2O2  which  was  usually 
<.  0.3  Torr. 

Total  error  in  the  determination  of  k^.  As  mentioned  earlier,  the  total 
estimated  error  in  calculating  the  fluence  is  12%.  If  we  add  to  this  error, 
possible  errors  in  the  knowledge  of  [H2O2],  the  quantum  yield  of  OH 
production,  and  the  uncertainties  in  absorption  cross  sections  of  H2O2  at  248.5 
and  202.6  nm  (5%  each),  the  total  uncertainty  in  our  calculation  of  [OH] 
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produced  upon  H2O2  photolysis  would  be  -15%.  It  is  important  to  note  here 
that  any  systematic  errors  in  the  knowledge  of  cr(H202,248.5  nm)  will  be 
canceled  out,  at  least  partially,  since  we  are  using  a(H2O2,202.6  nm)  and  as 
long  as  both  numbers  were  taken  from  the  same  set  of  data  the  ratio 
a(H202,248.5  nm)/(7(H2O2,202.6  nm)  will  be  much  more  accurate.  It  is  this 
ratio  that  is  important  in  calculating  [OHIq.  However,  we  have  conservatively 
assigned  5%  error  to  each  quantity. 

As  we  discussed  earUer,  not  all  OH  is  converted  to  HO2.  The  possible 
reduction  due  to  reaction  (21)  was  8%  in  the  worst  case.  In  addition,  the  HO2 
loss  due  to  reaction  (22)  varied  from  5%  to  9.7%  with  a  mean  value  of  7.5%. 
These  losses  have  been  well  accounted  for,  as  described  earlier.  However,  to 
be  on  a  generously  safe  side,  we  assume  that  the  errors  in  each  correction  could 
be  100%.  Using  all  these  errors  in  an  error  propagation  analysis  jdelds  the  net 
error  in  HO2  concentration  to  be  18%.  This  error  has  been  added  to  the 
precision  of  the  measurements  of  kj  and  the  total  estimated  errors  £u*e  shown 
in  the  last  column  of  Table  I. 
Comparison  with  Other  Studies 

Table  II  lists  all  measurements  of  kj  to  date.  This  table  is  essentially  a 
duplication  of  that  presented  by  Keyser^^  except  that  the  recent  results  of 
Sridharan  et  al.^^  and  those  of  the  present  investigation  are  included.  The 
possible  reasons  for  the  discrepancies  between  kj  values  measured  by  Keyser 
and  Sridharan  et  al.  on  one  hand  and  all  other  previous  measurements  on  the 
other  have  been  discussed  before^^'^^  and  need  not  be  repeated.  The  agreement 
between  our  results  and  that  of  Keyser,  (6.1  +.  0.4)  x  10"^^  cm^  molecule"^'^, 
is  extremely  good.  The  errors  quoted  above  for  Keyser's  value  do  not  include 
systematic  errors  which  he  estimates  to  be  +  25%.  The  results  of  Sridharan 
et  al.,^^  (5.4  ±  0.9)  x  10'^^  cm^  molecule'^s'^  are  also  in  good  agreement  with 
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our  results  considering  the  error  limits  of  the  two  studies.  (It  is  worth  pointing 
out  that  Keyser  and  Sridharan  et  al.  have  not  corrected  their  data  for  the 
occurrence  of  reaction  (22)  since  they  measured  [HO2]  and  the  change  in  0(^P) 
decay  rates  due  to  the  continuously  changing  concentration  of  HO2  will  be 
small  as  was  discussed  for  our  measurements.)  Based  on  the  k^  value  measured 
by  Keyser,  Sridharan  et  al.,  and  us,  it  appears  that  the  value  of  kx  is  well 
established  at  298  K  and  the  average  value  is  (5.9  +  1.0)  x  IQ-l^  cm^ 
molecule'^s"^,  independent  of  pressure.  The  atmospheric  implication  of  this 
higher  value  of  kj  compared  to  previous  results  have  been  discussed  by 
Keyser^^  and  Sridharan  et  al.^^  and  need  not  be  repeated. 
Comments  on  the  Pressure  Independence  of  k^ 

In  this  study  we  have  measured  ki  over  a  pressure  range  of  10  to  500  Torr 
of  N2,  a  factor  50  change.  A  plot  of  kj  as  a  function  of  pressure  is  shown  in 
Figure  10.  The  error  bars  that  are  shown  represent  precision  only.  The  reason 
for  doing  so  is  that  systematic  errors  are  essentially  the  same  at  all  pressures 
and  therefore  the  relative  values  of  kj  at  various  pressures  are  unaffected  by 
the  systematic  errors.  From  Figure  10,  it  is  clear  that  the  value  of  kj  is 
essentially  the  same  at  all  pressures  within  the  experimental  precision.  The 
maximum  range  of  kj  values  we  observe  is  5.81-6.47  x  10"^^  cm^  molecule"^s'^ 
The  average  value  falls  within  the  error  bounds  of  each  one  of  the  measured  k^ 
values  (including  that  in  100  Torr  of  Ar).  We  prefer  to  quote  a  pressure 
independent  value  of  (6.2  +  1.1)  x  10"^^  cm^  molecule-^s'l  for  kj.  Therefore,  it 
is  clear  and  worth  reiterating  that  the  value  of  kj  is  independent  of  pressure 
from  10  to  500  Torr. 

The  reaction  of  0(^P)  with  HO2  is  too  fast  to  be  understood  in  terms  of 
a  simple  H  atom  abstraction  reaction  proceeding  through  a  tight  transition 
state. ^^'^^  0(^P)  can  attack  HO2  either  on  the  H  atom  end  to  abstract  the  H 
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atom  or  on  the  oxygen  atom  end  to  form  an  H-0-0...0  (not  necessarily  lineeir) 
complex  which  can  decompose  to  give  OH  and  02^®.  The  latter  route  can 
proceed  on  a  purely  attractive  potential  energy  surface  since  the  free  radiced 
site  of  HO2  is  involved.  If  this  complex  is  formed,  it  might  be  possible  to 
stabiUze  it  and  decrease  its  decomposition  back  to  0(^P)  +  HO2.  If  this 
happens,  then  the  observed  value  of  k^  should  increase  at  high  pressures.  The 
absence  of  the  effect  of  pressure  on  kj  up  to  500  Torr  seen  here  indicates  that, 
if  formed,  the  lifetime  of  the  complex  towards  decomposition  into  0(^P)  +  HO2 
should  be  a  nanosecond  or  less.  This  short  lifetime  is  not  difficult  to 
comprehend  since  the  complex  would  be  very  energy  rich  and  the  energy  is 
distributed  over  relatively  few  degrees  of  freedom.  (A  discussion  of  subsequent 
experiments  that  clarify  the  reaction  pathway  is  presented  in  the  next  chapter.) 


CHAPTER  III 

TEMPERATURE  DEPENDENCE  OF  THE  0(3p)  +  HOg 
RATE  COEFFICIENT 


Introduction 
The  reaction  of  ground  state  oxygen  atoms  with  hydroperoxyi  radicals 

0(3p)  +  HO2  ^  OH  +  O2  (1) 

is  a  major  odd  oxygen  destruction  pathway  in  the  upper  stratosphere  and 
mesosphere.   Along  with  the  reactions 

0(3P)  +  OH  -  H  +  O2  (23) 

H  +  O2  +  M  ^  HO2  +  M  (25) 

H  +  O3  -*  OH  +  O2,  (26) 

reaction  (1)  plays  a  major  role  in  controlling  the  partitioning  among  H,  OH, 
and  HO2  radicals  in  the  upper  atmosphere.  Hence,  accurate  kinetic  data  for 
reaction  (1)  are  needed  in  order  to  model  upper  atmospheric  chemistry. 
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Several  kinetics  studies  of  reaction  (1)  are  reported  in  the  literature.^^'^®' 

23,36  Four  recent  direct  measurements  of  kj  at  298  K  are  in  excellent 
agreement,^^'^®'^^'^^  with  reported  rate  coefficients  all  within  the  range  (5.2-6.2) 
X  10'^^  cm^  molecule'^s'^.  However,  there  has  been  only  one  investigation  of 
the  temperature  dependence  of  kj.  Keyser^®  studied  reaction  (1)  over  the 
temperature  remge  229  -  372  K  in  a  discharge  flow  system  at  1  Torr  total 
pressure  and  observed  that  k^CT)  increased  with  decreasing  temperature;  his 
reported  "negative  activation  energy"  was  ~  0.4  kcal/mole.  Although  Keyser's 
study  of  reaction  (1)  was  a  high  quality  experiment  which  appears  to  be  free  of 
significant  systematic  errors,  the  uncertainty  in  his  reported  activation  energy 
will  remain  undesirably  high^^  until  independent  confirmation  is  reported. 

Several  years  ago,  we  developed  a  pulsed  laser  photolysis  technique  for 
studying  the  kinetics  of  radical-radical  reactions  at  pressures  up  to  1  atm.  We 
first  applied  this  technique  to  investigate  reaction  (1)  at  298  K  over  the 
pressure  range  10  -  500  Torr.^^  In  this  paper,  we  report  the  results  of  a  pulsed 
laser  photolysis  resonance  fluorescence  study  of  the  temperature  dependence 
of  kj.  Our  results,  obtained  using  a  much  different  experimental  approach  than 
that  employed  by  Keyser,^^  confirm  his  reported  temperature  dependence  and 
also  demonstrate  that  k^  is  independent  of  pressure  at  subambient 
temperatures. 

Experimental  Section 

With  a  few  modifications,  the  experiments  were  carried  out  in  the  same 
manner    as   our   previous   room   temperature   study. ^^      A   review   of  the 
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experimental  approach,  along  with  details  pertinent  to  this  investigation,  is 

given  below.   A  schematic  of  the  apparatus  is  shown  in  Figure  11. 

All  experiments  were  carried  out  under  "slow  flow"  conditions  using  a 
jacketed  Pyrex  reactor  with  an  interned  volume  of  320  cm^.  The  cell  was 
medntained  at  a  constant  temperature  by  circulating  ethylene  glycol  from  a 
thermostated  bath  through  the  outer  jacket.  A  copper-constantan  thermocouple 
with  a  stainless  steel  jacket  was  inserted  into  the  reaction  zone  through  a 
vacuum  seal,  thus  allowing  measurement  of  the  gas  temperature  under  the 
precise  pressure  and  flow  conditions  of  the  experiment. 

Reactants  were  produced  with  HO2  in  excess  using  248.5  nm  pulsed  laser 
photolysis  of  H2O2/O3/N2  mixtures: 


H2O2  +  hv(248.5nm)  ^  2  OH 


(13) 


O3  +  hv(248.5nm)- 


^  0(lD)  +  OaCalAJ 


^  0(3P)  +  OgCX^S -) 


(14a) 


(14b) 


OH  +  H2O2  -*  HO2  +  H2O 


(15) 


0(lD)  +  N2  -  0(3P)  +  N2 


(16) 
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A  Lambda  Physik  model  200E  KrF  excimer  laser  was  used  as  the  photolysis 

light  source.  Kinetic  data  was  obtained  by  time-resolved  resonance  fluorescence 
detection  of  0(^P).  The  laser  pulsewidth  was  20  ns  while,  under  our 
experimental  conditions,  reactions  (15)  and  (16)  proceeded  at  rates  of  (3-10)  x 
10^  and  7  x  lO'^s'^,  respectively.  0(^P)  decay  rates  were  typically  in  the  remge 
30  -  500  s-1. 

As  in  our  previous  study,  the  concentration  of  the  excess  reactant,  HO2, 
was  not  directly  measured  but  was  calculated  based  on  experimentally 
measured  and  other  known  parameters.  To  obtain  [HO2]  one  must  determine 
[OH]o,  the  initial  concentration  of  photolytically  produced  OH,  and  the  yield  of 
HO2  from  reaction  (15)  and  competing  side  reactions.  The  chemistry  of 
conversion  of  OH  to  HO2  is  discussed  in  detail  in  a  later  section.  Since  the 
concentrations  of  H2O2  and  O3  were  such  that  the  system  was  optically  thin  at 
248.5  nm,  [OHIq  could  be  calculated  from  the  following  relationship 

[OHJo  =  <DoH  X  a(H202,248.5  nm,T)  x  [H2O2]  x  F  (X) 

where  <I>oh  is  the  quantum  yield  for  OH  production  from  248.5  nm  photolysis 
of  H2O2,  a(H202,248.5  nm,T)  is  the  absorption  cross  section  for  H2O2  at  248.5 
nm  and  temperature  T,  and  F  is  the  laser  photon  fluence.   The  determination 
of  each  factor  in  equation  (I)  is  discussed  in  detail  below. 
^OH-   It  is  known25>26,38  ^hat  ^qh  =2. 
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a(H202-248.5  nm.T).  The  absorption  cross  section  for  H2O2  at  248.5  nm 

is  known  to  be  8.8  x  lO'^^  cm^  at  298  K.^*^     We  have  recently  measured 

temperature  dependent  absorption  cross  sections  for  hydrogen  peroxide  over 

the  wavelength  range  193-350  nm.^^    At  248.5  nm,  the  following  expression 

describes  the  observed  temperature  dependence  (units  are  cm^  molecule"^): 

a(H202,248.5  nm,  T)  =  1.023  x  lO'^^  exp{(-45  ±  20)/T}  (XI) 

IH2Q21-  Hydrogen  peroxide  can  be  lost  in  the  slow  flow  system  either  by 
decomposition  (particularly  at  higher  temperatures)  or  by  condensation  (at 
lower  temperatures).  To  ensure  that  the  H2O2  concentration  in  the  reactor  was 
known,  we  monitored  H2O2  by  UV  photometry  before  the  gas  mixture  entered 
the  reactor  and  after  the  gas  mixture  exited  the  reactor.  The  absorption  cells 
were  216.2  and  90.0  cm  in  length.  The  monitoring  wavelengths  were  228.8  nm 
in  the  longer  cell  (Cd  line)  and  202.6  nm  in  the  shorter  cell  (Zn"*"  line).  Both 
cells  were  kept  at  ambient  temperature,  and  the  absorption  cross  sections 
needed  to  convert  absorbance  data  to  H2O2  concentration  were  obtained  by 
interpolation  from  current  NASA  recommendations:^*^  1.86  x  10"^®  cm^  at  228.8 
nm  and  4.31  x  10'^^  cm^  at  202.6  nm.  When  the  reactor  temperature  was  350 
K  or  below,  the  difference  in  H2O2  concentration  measured  in  the  two 
absorption  cells  was  never  more  than  a  few  percent.  At  391  K,  the  highest 
temperature  at  which  experiments  were  performed,  15-20%  of  the  H2O2  was 
lost  upon  traversal  of  the  reactor.  The  H2O2  concentration  in  the  reaction  zone 
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was  always  taken  to  be  the  temperature-corrected  average  of  the  concentrations 

measured  in  the  two  absorption  cells.  Under  our  experimental  conditions, 
absorption  by  ozone  was  negligible  compared  to  absorption  by  H2O2  at  202.6 
nm.  At  228.8  nm,  ozone  made  a  small  but  significant  contribution  to  the  total 
absorbance.  Care  was  taken  to  ensure  that  [O3]  was  constant  during  the  I  and 
Iq  measurements  required  for  the  [H2O2]  determination. 

F.  The  photolysis  laser  beam  was  made  spatially  uniform  through  use  of 
a  segmented  aperture  optical  integrator.  ^^'^^  Virtually  the  whole  cross-sectional 
area  of  the  cell  was  irradiated,  and  the  depth  of  focus  of  the  integrated  beam 
was  such  that  radical  concentrations  were  nearly  uniform  down  the  entire 
length  of  the  cell.  The  laser  beam  fluence  was  measured  as  the  besun  exited  the 
reactor  by  using  an  EG&G  photodiode  based  radiometer  capable  of  measuring 
individual  pulses.  Pulse-to-pulse  stability,  a  requirement  for  signal  averaging 
in  this  type  of  experiment,  was  found  to  be  very  good.  Only  an  occasional  pulse 
energy  deviated  from  the  average  by  more  than  ±_  5%.  In  order  to  avoid  having 
to  correct  the  measured  fluence  for  reflection  off  the  back  window,  an 
antireflection  coated  window  with  >  99.5%  transmission  at  248.5  nm  was 
employed.  The  radiometer  was  calibrated  by  using  a  novel  ozone  actinometry 
method  which  has  been  described  in  detail  previously.  ^^ 

As  mentioned  above,  all  experiments  were  carried  out  under  slow  flow 
conditions.  The  linear  flow  velocity  through  the  reactor  was  typically  12  cm  s'^ 
and  the  laser  repetition  rate  was  typically  0.4  Hz.  The  reactor  was  23  cm  in 
length,  so  the  reactor  volume  was  completely  replenished  with  a  fresh  gas 
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mixture  between  laser  pulses.  All  experiments  employed  nitrogen  as  the  buffer 

gas  at  a  total  pressure  of  80  Torr.  Data  were  obtained  over  the  temperature 
range  266-391  K.  The  temperature  range  was  limited  at  the  low  end  by  the 
vapor  pressure  of  H2O2  and  at  the  high  end  by  the  thermal  instability  of  H2O2. 

The  gases  used  in  this  study  had  the  following  stated  minimum  purities: 
N2,  99.999%;  O2,  99.99%.  Hydrogen  peroxide  was  90  wt  %  in  water.  It  was 
concentrated  further  by  bubbling  N2  through  the  semiple  for  several  days  before 
experiments  were  undertaken  and  continuously  during  the  course  of  the 
experiments.  To  prevent  significant  decomposition  of  H2O2,  all  components 
traversed  by  H2O2  between  the  bubbler  and  the  exit  from  the  last  absorption 
cell  were  Pyrex  or  Teflon  with  the  exception  of  a  few  stainless  steel  fittings. 
The  needle  valve  and  flowmeter  in  the  H2O2  line  were  positioned  so  that  N2 
flowed  through  these  components  before  entering  the  bubbler.  Ozone  was 
prepared  by  passing  O2  through  a  commercial  ozonator  and  was  stored  on  silica 
gel  at  195  K.  Before  use  it  was  degased  at  77  K  to  remove  O2.  Dilute  O3/N2 
mixtures  were  prepared  in  12  liter  Pyrex  bulbs  for  use  in  experiments. 

A  typical  experiment  was  initiated  by  flowing  (3-10)  x  10^^  O3  molecules 
cm"^  in  80  Torr  of  N2  through  the  reactor  and  absorption  cells.  Iq  and  J^,  the 
intensities  of  228.8-  and  202.6-nm  light  transmitted  through  the  absorption 
cells,  were  measured.  Next,  H2O2  was  introduced  into  the  gas  flow,  the  N2  flow 
was  reduced  so  the  total  flow  rate  and  total  pressure  were  the  same  after 
addition  of  H2O2  as  before  addition  of  H2O2,  and  I  and  J,  the  reduced 
intensities  of  228.8-  and  202.6  nm  light  transmitted  through  the  reactor,  were 
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measured.   Using  the  measured  values  of  I,  Iq,  J,  and  J^,  the  concentration  of 

H2O2  in  the  gas  stream  entering  and  exiting  the  reactor  was  calculated;  it 
ranged  from  (2-6)  x  10^^  molecules  cm'^.  I  and  J  were  continuously  monitored 
during  the  course  of  an  experiment.  The  multichannel  analyzer  was 
pretriggered  before  the  photolysis  laser  fired  to  obtain  the  background  count 
rate.  The  concentration  of  0(^P)  was  monitored  as  a  function  of  time  after  the 
photolysis  pulse.  A  total  of  50  -  200  laser  shots  were  averaged  to  obtain  one 
pseudo-first  order  kinetic  decay.  The  fluence  of  each  laser  pulse  was  measured 
using  the  radiometer  and  a  calibrated  aperture.  The  laser  fluence  was  varied 
at  constemt  [O3]  and  [H2O2]  to  obtain  pseudo-first  order  decays  as  a  function 
of  [HO2].  Five  to  ten  fluence  values  were  employed  to  determine  each 
bimolecular  rate  constant,  k^,  from  the  slope  of  a  k'  versus  [HO2]  plot  (k'  =  the 
0(^P)  pseudo-first  order  decay  rate).  After  the  fluence  variations  were 
completed,  the  H2O2  flow  was  turned  off,  the  total  pressure  and  total  flow  rate 
were  readjusted,  and  Iq  and  J^,  were  remeasured.  All  HO2  concentrations 
employed  in  the  k^  determinations  were  in  the  range  (0.15-8.5)  x  10^^  molecules 
cm'^.  It  should  be  noted  that,  in  a  series  of  runs  involving  variation  of  the  laser 
fluence  at  constant  [O3]  and  [H2O2],  the  ratio  [OH]o/[0(^P)]o  is  not  altered;  to 
change  this  ratio,  the  composition  of  the  reaction  mixture  must  be  varied. 

Results  and  Discussion 
In  the   absence   of  competing  side   reactions   which   affect   the   HO2 
concentration,  the  kinetic  system  is  inherently  pseudo-first  order,  i.e.,  all  HO2 
lost  via  reaction  (1)  is  rapidly  regenerated  via  reaction  (15).   However,  as  will 
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be  discussed  in  detail  below,  the  importance  of  certain  side  reactions  is 

suppressed  when  [HO2]  >  >  [0(^P)].  Most  of  our  experiments  were  carried  out 

under  experimental  conditions  where  [HO2]  «  10[O(^P)lo,  although  this  ratio 

was  varied  as  a  check  on  the  kinetic  model  xised  to  extract  k^.     At  low 

temperature,  where  relatively  low  H2O2  concentrations  had  to  be  employed,  the 

HO2  to  0(^P)  ratio  was  typically  somewhat  lower  than  that  employed  at  higher 

temperatures. 

Photolytically  produced  0(^P)  can  be  lost  via  the  following  pseudo-fu-st 

order  processes: 

0(3P)  +  HO2  ^  OH  +  O2  (1) 

0(3P)  +  H2O2  ^  OH  +  HO2  (17) 

0(3P)  +  03^  2O2  (6) 

0(3p)  ^  loss  by  diffusion  from  the  detector  field  of 

view  and  reaction  with  background  impurities  (18) 

Under  our  experimental  conditions,  reaction  (1)  dominated  0(^P)  removal 
exccept  at  very  low  HO2  levels.  Reaction  (6)  was  of  negligible  importance  while 
reaction  (17)  was  minor  but  not  negligible. 3^.    The  value  of  kjg  was  directly 
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measured  to  be   -5  s '^  1  to  2  orders  of  magnitude  slower  than  ki[H02]  under 

most  experimental  conditions. 

In  the  absence  of  competing  side  reactions  which  affect  the  concentrations 
of  0(^P)  or  HO2,  removal  of  0(^P)  should  obey  first  order  kinetics 

ln{[0(3p)]o/[0(3p)]}  =  (kiCHOg]  +  kd)t  =  k'e^ptit  (XII) 

where 

kd  =  knLHaOa]  +  keLOg]  +  kig  (XIII) 

Figure  12  shows  typical  plots  of  In  [0(^P)]  versus  time.  Equation  XII  does 
indeed  appear  to  be  obeyed.  If  side  reactions  were  unimportant,  the 
bimolecular  rate  coefficient  would  be  obtained  from  the  slope  of  a  k'^^ptj  versus 
[OH]o  plot,  such  as  shown  in  Figure  13. 

Secondary  Chemistry 

In  the  paper  describing  our  298K  study  of  reaction  (1),  possible  side 
reactions  were  considered  in  detail.^^  It  was  concluded  that  two  reactions  could 
significantly  affect  the  HO2  temporal  profile: 

OH  +  HO2  -*  O2  +  H2O  (21) 

HO2  +  HO2  ->  H2O2  +  O2  (22) 
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Figure  12. 


Typical  0(^P)  temporal  profiles  observed  following  248.5 
nm  pulsed  laser  photolysis  of  O3/H2O2/N2  mixtures. 
Experimental  conditions:  T  =  281  K,  P  =  80  Torr, 
[H2O2]  =  3.9  X  IQiS  molecules  cm-3,  [O3]  =  6  x  10^2 
molecules  cm'^,  laser  fluence  (in  units  of  mJ  cm'^)  =  (a) 
1.95,  (b)  4.55,  and  (c)  6.31.  Solid  lines  are  obtained 
from  least  squares  analyses  of  the  first  two  1/e  times  of 
0(^P)  decay  and  give  the  following  pseudo-first  order 
decay  rates  (k'^^pti):  (a)  HI  s'l,  (b)  210  s'l,  (c)  303  s"!. 
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Figure  13.  Typical  plots  of  k'gjpti  versus  [OHIq  (open  circles)  and  k' 
versus  [H02]inas  (dosed  circles).  Experimental  conditions:  T 
=  281  K,  P  =  80  Terr,  [HgOg]  =  3.9  x  lO^^  molecules  cm-^, 
[O3]  =  6  X  10^^  molecules  cm'"^.  The  dashed  line  is  obtained 
from  a  linear  least  squares  analysis  of  the  k'^^-ti  versus  [OHlg 
data  and  gives  the  "uncorrected"  rate  coefficient  (5.24  ±_  0.38) 
X  10"^^  cm^  molecule"^  s'^.  The  solid  line  is  obtained  from  a 
linear  least  squares  analysis  of  the  k'  versus  [H02]inax  ^^^^ 
and  gives  the  "corrected"  rate  coefficient  (6.15  jf  0.54)  x  10"^^ 
cm^  molecule'^  s'^. 
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To    quantify   the   roles   of  reactions    (21)    and    (22)    in   our   kinetics 

experiments,  a  series  of  computer  simulations  were  carried  out  where  the 
temporal  profiles  of  key  species  were  calculated  under  a  variety  of  experimental 
conditions  by  numerical  integration  of  the  appropriate  rate  equations.  (The 
final  section  of  Appendix  A  contains  a  discussion  of  the  integration  procedure.) 
For  completeness,  a  number  of  reactions  which  were  expected  to  play  very 
minor  roles  in  the  0(^P)  and  HO2  kinetics  were  included  in  the  mechanism. 
The  complete  set  of  reactions  and  rate  coefficients  used  in  the  simulations  is 
given  in  Table  III. 

Reaction  (21)  competes  with  reaction  (15)  during  the  period  immediately 
after  the  laser  pulse  when  OH  is  being  converted  to  HO2.  Each  time  reaction 
(21)  occurs,  two  HO2  radicals  are  lost  which  otherwise  would  have  been  present 
to  react  with  0(^P).  Reaction  (21)  is  most  important  at  high  laser  fluence  (i.e., 
high  [OH]o/[H202])  and  at  low  temperature.  From  the  computer  simulations, 
a  set  of  curves  were  constructed  which  relate  [H02]njax>  ^^^  peak  amount  of 
HO2  present  after  all  photoljiiically  produced  OH  has  reacted  away  but  before 
appreciable  HO2  loss  via  processes  such  as  reaction  (22)  has  occurred,  to  [OHIq. 
Representative  correction  curves  are  plotted  in  Figure  14.  For  the  101 
experiments  used  in  the  ki(T)  determinations,  the  average  value  for 
[H02]niax/[C)H]o  was  0.947  while  the  minimum  value  under  any  set  of 
experimental  conditions  was  0.845. 

As  pointed  out  above,  loss  of  0(^P)  via  reactions  (1),  (6),  (17)  and  (18) 
occurs  on  a  time  scale  which  is  long  compared  to  the  time  scale  for  HO2 
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Table  III.   Reaction  Set  for  Computer  Simulations. 


Reaction 

Rate  Coefficient<^^) 

1          O  +  HO2  -  OH  +  O2 

ki  =  3  X  10-11  exp(200/T) 

0  +  OH  -  H  +  O2 

k23  =  2.2  X  10-11  exp(117/T) 

H  +  O3  -  OH  +  O2 

k26  =  1.4  X  10-10  exp(-470/T) 

OH  +  H2O2  -*  HO2  4-  HgO 

ki5  =  3.1  X  10-12  exp(-187/T) 

0  +  H2O2  ^  OH  +  HO2 

ki7  =  1.4  x  10-12  exp(-2000/T) 

0  +  O3  ^  2  O2 

kc  =  8  x  10-12  exp  (-2060/T) 

0  -♦  loss 

ki8  =  5  s-l 

OH  +  HO2  -  O2  +  H2O 

k2i  =  1.7  X  10-11  exp(416/T)  +  3  x 
10-31  [M]  exp(500/T) 

HO2  +  HO2  -  H2O2  +  O2 

k22  =  2.3  X  10-13  exp(590/T)  +  1.7  x 
10-33[M]  exp(1000/T) 

HO2  -*■  loss 

k27  =  5  s-l 

HO2  +  O3  -  OH  +  2  O2 

kg  =  1.4  X  10-14  exp(-580/T) 

H  +  HO2  -►  2  OH 

k24  =  6.4  X  10-11  exp(0/T)              1 

(a)       All  rate  coefficients  are  taken  from  reference  37  except  k^s  which  was 
measured  and  k27  which  was  set  equal  to  kjs- 

Ob)       All  rate  coefficients  except  k^g  and  k27  are  in  units  of  cm^  molecule-1  s-l. 
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1.00- 


0.80 


10 


[OH]o  (lO^^molecules  per  cm^) 


Figure  14.  Typical  correction  curves  (obtained  from  computer 
simulations)  which  relate  [H02]inax  <^o  [OHIq.  All  curves 
shown  in  the  figure  are  from  simulations  with  [H2O2]  = 
4  X  1015  molecules  cm'^.  Temperature:  I,  400  K;  II,  300 
K;  III,  260  K.    [OH]o/[0(3p)]o:  a,  lO"*;  b,  10;  c,  4. 
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formation.    Hence,  if  the  concentration  of  HO2  were  constant  over  the  time 

period  of  0(^P)  removal,  kj  could  be  obtained  from  the  slope  of  a  k'^^^i  versus 
[H02]niax  plot.  Unfortunately,  small  but  significant  time  variation  of  [HO2]  can 
result  from  the  occurrence  of  reactions  (17),  (22),  and  (27).  Reactions  (17)  and 

HO2  -*  loss  by  diffusion  from  the  reaction  zone  and 

reaction  with  background  impurities  (27) 

(27)  are  most  important  at  low  laser  fluences,  i.e.,  low  [HO2].  Reaction  (17)  is 
strongly  temperature  dependent  and  increases  in  importance  at  high 
temperature.  Reaction  (17)  is  also  a  more  important  source  of  HO2  when 
[OH]o/[0(^P)]o  is  relatively  low.  Reaction  (22)  becomes  an  important  HO2 
removal  mechanism  at  relatively  high  HO2  concentrations.  To  avoid  large 
corrections  for  HO2  loss  via  reaction  (22),  all  experiments  were  carried  out  with 
[H02],nax  <  9  X  10^^  molecules  cm"^  and  all  data  analyses  were  restricted  to 
two  1/e  times  of  0(3p)  decay;  i.e.,  no  data  where  [0(3P)]/[0(3P)]q  <  0.13  were 
used  in  the  data  analysis.  Computer  simulations  were  carried  out  under  a 
variety  of  experimental  conditions.  The  first  two  1/e  times  of  the  slightly 
nonexponential  computer  generated  0(^P)  temporal  profiles  were  least  squares 
fit  to  an  exponential  decay  to  obtain  k'gj^^,  the  simulated  decay  rate.  The 
simulated  decay  rates  were  then  compared  to  the  "real"  decay  rates,  k'^,  to 
obtain  a  set  of  correction  curves,  some  of  which  are  shown  in  Figure  15. 
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k',  =  kiLHOg]^  +  knCHgOa]  +  kgLOg]  +  kig  (XIV) 

The  first  two  1/e  times  of  each  experimental  0(^P)  temporal  profile  were  least 
squares  fit  to  equation  (XII)  to  obtain  k'^jp^j.  A  corrected  value  of  k'  was  then 
computed  from  the  expression 

k'  =  k'e^ptl(k'r/k'8im)  (XV) 

One  source  of  uncertainty  in  the  computer  simulations  concerns  the  fact  that 
k27  was  not  measured;  it  was  estimated  that  k27  =  k^g  ~  5  s'^.  Lack  of 
knowledge  of  the  exact  value  of  k27  increases  the  uncertainty  in  the  k'r/k'gjj^ 
factors  at  low  concentrations  of  HO2.  The  low  [HO2]  data  are  relatively 
unimportant  in  defining  kj,  so  the  uncertainty  of  k27  makes  only  a  minor 
contribution  to  the  overall  undertainty  in  kj.  For  the  101  experiments  used  in 
the  ki(T)  determinations,  the  average  value  of  ^'r/^'eim  ^^^  ^-^^  ^^^  ^^^ 
maximum  value  under  any  set  of  experimental  conditions  was  1.18. 

As  discussed  above,  k^  values  corrected  for  secondary  chemistry  were 
obtained  from  the  slopes  of  plots  of  k'  versus  [H02]niax-  TVpical  data  are 
shown  in  Figure  13.  For  all  15  rate  coefficients  measured,  the  corrected  k^ 
values  were  larger  than  the  values  which  would  have  been  obtained  from  plots 
of  k'e^ti  versus  [OHlg  (see  Figure  13,  for  example).  The  magnitude  of  the 
secondary  chemistry  corrections  was  largest  at  the  lowest  temperatures.  For 
the  15  rate  coefficients  reported,  the  average  difference  between  the  corrected 
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and  uncorrected  bimolecular  rate  coefficients  was   16%  while  the  largest 

difference  was  26%. 
Summary  of  Results 

The  experimental  results  are  summarized  in  Table  IV.  Errors  quoted  for 
individual  "corrected"  k^  determinations  are  2a  and  refer  only  to  the  precision 
of  the  k'  versus  [H02]niax  data.  The  absolute  accuracy  of  our  k^  determinations 
is  limited  not  only  by  precision  but  also  by  uncertainties  in  measurement  of  the 
laser  photon  fluences  (F),  the  H2O2  concentration,  the  correction  factor  used 
to  obtain  [H02]inax  from  [OHIq,  the  correction  factor  used  to  obtain  k'  from 
k'expti  (i.e.,  k'r/k'gijjj),  and  unidentified  systematic  errors.  We  estimate  the 
pertinent  2a  uncertainties  to  be  as  follows:  F,  10%;  [H2O2],  5%;  [H02]niax/[OHlo 
5%;  k'y/k'gim.  10%  at  266K  and  5%  at  298K  and  above;  unidentified  systematic 
errors,  5%.  Precision  did  not  appear  to  be  temperature  dependent  and 
averaged  9%.  Hence,  the  absolute  accuracy  of  an  individual  kj  determination 
is  estimated  to  be  +  23%  at  266  K  and  +  21%  at  298  K  and  above. 

An  Arrhenius  plot  of  our  results  is  shown  in  Figure  16.  An  unweighted 
linear  least  squares  analysis  of  the  In  k^  versus  1/T  data  gives  the  following 
expression  (units  are  cm^  molecule"^  s'^): 

ki(T)  =  (2.91  ±  0.70)  X  10-11  exp[(228  ±  75)/T]  (XVI) 
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Figure  16.     Arrhenius  plot  for  the  0(^P)  +  HO2  reaction. 
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The  uncertainties  quoted  in  the  above  expression  are  2a  and  represent  precision 

only. 

Comparison  with  Previous  Work 

Our  results  are  compared  with  those  reported  by  other  investigators  in 
Table  V.  The  298  K  value  reported  in  this  paper  is  identical  with  our 
previously  reported  value  ^^  and  at  the  upper  end  of  the  group  of  recent  direct 
determinations^^'^^'^^'^®  which  span  the  range  (5.2-6.2)  x  10'^^  cm^  molecule'^ 
s'^.  Other  than  our  studies,  the  other  recent  direct  measurements  all  employed 
low  pressure  discharge  flow  systems,  all  measured  not  only  0(^P)  but  eJso  HO2 
(either  directly  or  indirectly  by  conversion  to  OH),  and  all  carefully  considered 
the  role  of  competing  side  reactions.  Hence,  there  is  no  obvious  reason  to 
prefer  one  value  over  another.  The  early  work  of  Burrows  et  al.^^  and  Hack 
et  al.^^  give  ki(298  K)  values  which  are  considerably  lower  than  those  reported 
in  references  16,  19,  20,  and  36.  Possible  reasons  for  the  apparently  erroneous 
results  reported  in  references  22  and  21  are  discussed  elsewhere. ^^  Lii  et  al.^^ 
obtained  a  value  for  k^  by  comparing  computer  simulations  with  HO2  and  O3 
concentration  profiles  observed  following  pulsed  radiolysis  of  02/H2/Ar 
mixtures.  Their  reported  rate  coefficient  agrees  well  with  our  results. 
However,  as  pointed  out  by  Keyser,^^  the  experimental  data  of  Lii  et  al.  are 
very  insensitive  to  the  value  of  kj;  hence,  their  error  limits  should  be  ~500% 
rather  than  the  reported  30%. 

Very  little  temperature  dependent  data  are  available  for  reaction  (1). 
Values  for  kj  of  8  x  lO'll  cm^olecule-^s-l  at  1600  K^^  and  6  x  lO'l^cm^ 
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molecule'^s"^  at  1050  K'^^  have  been  inferred  from  flame  studies.  The  only 
previous  temperature  dependence  study  in  the  atmospheric  temperature  regime 
is  that  of  Keyser.^^  As  seen  from  the  comparisons  in  Table  V  and  Figure  14, 
our  results  are  in  excellent  agreement  with  those  of  Keyser.  Since  our 
experiments  were  very  different  in  methodology  from  Keyser's  and  are  subject 
to  different  sources  of  systematic  errors,  the  uncertainty  of  k^CT)  for 
atmospheric  modeling  purposes  is  now  greatly  reduced.  Keyser's  study 
employed  helium  buffer  gas  at  a  pressure  of  1.0  Torr  while  our  study  employed 
nitrogen  buffer  gas  at  a  pressure  of  80  Torr.  Thus,  agreement  between  the  two 
studies  strongly  suggests  that  kj  is  independent  of  pressure  over  the  relevant 
upper  atmospheric  temperature  and  pressure  ranges. 

Reaction  (1)  could  proceed  via  a  hydrogen  abstraction  mechanism  (la)  or 
via  formation  of  an  energized  HOOO  "pseudo-intermediate"  (lb): 

O  +  HO2  -*  O HO2  ^  OH  +  O2  (la) 

O  +  O2H  -*  O    OOH  -*  O2  +  OH  (lb) 

Sridharan  et  al.'^^  recently  reported  an  elegant  experiment  in  which  ^^OH  and 
^®0H  products  from  the  reaction  of  ^^0  with  HO2  were  monitored  in  a 
discharge  flow  system.  They  found  that  only  ^^OH  was  produced,  implying  that 
O  reacts  with  HO2  via  channel  (lb).  Thermochemical  estimates  suggest  that 
HOOO  is  bound  relative  to  O  +  HO2,  but  is  ""  12  kcal  mol'^  less  stable  thein  the 
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OH  +  O2  products.'^^''^  Mozurkewich'*^  points  out  that  since  HOOO  is  bound 
relative  to  O  +  HO2,  we  might  expect  to  find  a  long-range  interaction  that 
would  produce  a  transition  state  for  reaction  (1)  very  similar  to  that  expected 
from  formation  of  HOOO.  Mozurkewich's  RRKM  calculations  yield  a  rate 
constant  for  formation  of  HOOO  of  5.5  x  10"^^cm^olecule'^s"^  independent  of 
temperature.'^^  The  observed  negative  activation  energy  for  reaction  (1)  could 
probably  be  reproduced  if  a  small  barrier  were  assumed  in  the  HOOO  -*  OH 
+  OH2  reaction  path.'*^ 
Implications  for  Atmospheric  Chemistry 

Model  calculations  of  OH  and  HO2  concentration  profiles  in  the  upper 
stratosphere  are  very  sensitive  to  the  choice  of  k^CT).  Kaye  and  Jackman,^*^ 
considering  both  sensitivity  of  their  model  to  variovis  parameters  and 
uncertainties  in  these  parameters,  have  concluded  that  the  uncertainty  of  kj 
contributes  more  to  the  uncertainty  in  [OH]  and  [HO2]  at  35°N,  40  km  altitude 
than  any  parameter  in  their  model  except  k2i.  The  currently  recommended 
value37  for  kjCT)  is 

ki(T)  =  3.0  x  10-llexp[(200+200)/T]  cm^molecule-ls"!  (XVII) 

The  results  reported  in  this  paper  will  have  little  effect  on  the  recommended 
A  factor  and  activation  energy,  but  will  substantially  reduce  the  uncertainty  in 
the  above  expression;  this  vnll,  in  turn,  significantly  reduce  the  overall  uncer- 
tainty in  model  calculations  of  upper  stratospheric  OH  and  HO2  concentrations 
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and  facilitate  meaningful  comparisons  of  stratospheric  measurements  with 

photochemical  models. 

Jackman  et  al.'*®  have  recently  compared  O3  concentrations  at  5°N  and  43 
km  altitude  measured  by  LIMS  (limb  infrared  monitor  of  the  stratosphere)  with 
those  calculated  from  a  photochemical  model  using  LIMS  measurements  of 
H2O,  HNO3,  NO2,  and  temperature  and  SAMS  (stratospheric  and  mesospheric 
sounder)  measurements  of  CH4  as  input.  The  model  predicted  lower  O3  levels 
than  those  actually  observed.  A  sensitivity  analysis  showed  that  the  overall 
uncertainty  in  the  model  calculation  was  a  factor  of  1.7,  about  the  same 
magnitude  as  the  discrepancy  between  model  and  measurement.  Of  the  many 
parameters  used  in  the  model,  the  uncertainty  of  kj  was  found  to  make  the 
sixth  largest  contribution  to  the  overall  uncertainty  of  the  calculation.  Our 
results  will  therefore  result  in  a  small  but  significant  reduction  in  the  model 
uncertainty. 


CHAPTER  IV 

PULSED  LASER  PHOTOLYSIS  KINETICS  STUDY 
OF  THE  0(3P)  +  CIO  REACTION 


Introduction 
The  reaction  of  ground  state  oxygen  atoms  0(^P)  with  CIO  radicals  is  the 
rate  determining  step  in  the  dominant  catalytic  cycle  via  which  chlorine  atoms 
destroy  odd  oxygen  in  the  middle  stratosphere: 

O  +  CIO  -*  O2  +  CI  (2) 

CI  +  O3  -  CIO  +  O2  (12) 


O  +  O3  -  2O2  net 


The   primary   source   of  stratospheric   chlorine   atoms   is   the   photolysis   of 
anthropogenic  chlorofluorocarbons. 

Seven  measurements  of  k2(298  K)  are  reported  in  the  literature."*®"^^  There 
is  agreement  among  the  five  most  recent  studies  that  k2(298  K)  lies  in  the  range 
3.5  -  4.2  X  10'^^  cm'^olecule'^s"^.  The  activation  energy  for  reaction  (2)  is 
known  to  be  small,^^"^^  but  its  value  is  not  as  well  defined  as  would  be  desirable 
for  such  an  important  stratospheric  reaction.    In  fact,  it  is  not  clear  if  k2(T) 
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increases   or   decreases   with  decreasing  temperature.      In   addition  to   the 

abovementioned  studies  of  reaction  (2)  at  atmospheric  temperatures,  one  high 

temperature  (1250  K)  shock  tube  measurement  of  kj  has  been  reported,^*^  as  has 

one  theoretical  calculation  of  k2(T)  over  the  temperature  range  220-1000  K  ^ 

Because  predictions  of  chlorine  catalyzed  ozone  loss  are  very  sensitive  to 
the  value  of  k2(T)  used  in  model  calculations,  it  is  important  that  this  rate 
coefficient  be  determined  with  high  precision  at  stratospheric  temperatures. 
Studies  employing  a  variety  of  experimental  techniques  are  desirable  in  order  to 
uncover  possible  systematic  errors.  All  previous  studies  of  reaction  (2)  at 
ambient  and  subambient  temperatures^^"^^  employed  discharge  flow  systems 
which  were  limited  to  toted  pressures  of  10  Torr  or  less.  It  is  interesting  to  note 
that  reaction  (2)  occurs  on  a  potential  energy  surface  with  a  minimum  along  the 
reaction  coordinate,  i.e.,  the  intermediate  complex  ClOO  is  a  bound  species 
whose  ground  state  correlates  with  0(^P)  +  CIOCX^jt)^^.  Reactions  which  occur 
on  potential  energy  surfaces  of  this  type  often  exhibit  negative  activation 
energies  and  pressure  dependent  rates. 

We  have  recently  developed  a  pulsed  laser  photolysis  method  for  carrying 
out  direct  kinetics  studies  of  radical-radical  reactions  at  pressures  up  to  one 
atmosphere,  and  applied  this  method  to  study  the  temperature  and  pressure 
dependences  of  the  O  +  HO2  reaction. ^^'^"^  Using  an  extension  of  the  technique 
employed  in  the  O  +  HO2  investigations,  we  have  studied  the  kinetics  of  reaction 
(2)  in  N2  buffer  gas  over  the  temperature  and  total  pressure  ranges  231-367  K 
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and  25-500  Torr.  Our  results,  which  include  observation  of  a  significant  negative 

activation  energy,  are  reported  in  this  paper. 

Experimental 

A  schematic  of  the  apparatus  appears  in  Figure  17.  The  two  radical  species 

were  created  via  a  scheme  involving  two  separate  photolysis  lasers.  Under  "slow 

flow"  conditions,  a  gas  mixture  containing  CI2  and  O3  in  a  large  excess  of  N2 

buffer  gas  was  first  subjected  to  photolysis  by  a  XeF  excimer  laser, 

CI2  +  hv   (351  nm)  -*  2C1  (28) 

The  combination  of  the  excimer  laser  fluence  and  [CI2]  was  always  large  enough 
for  the  condition  [CUq  >  [03]o  to  hold.  During  a  predetermined  delay  period, 
t',  the  reaction 

CI  +  O3  -►  CIO  +  O2  (12) 

was  allowed  to  go  to  completion.  At  this  time  the  ozone  in  the  reaction  cell  had 
effectively  been  titrated  by  CI  atoms  and  the  initial  value  of  [03]^,  could  be 
related  to  [ClO]^'.  At  the  end  of  this  delay  a  second  laser  pulse,  the  fourth 
harmonic  of  the  fundamented  wavelength  from  a  Nd:YAG  laser,  photolyzed  a 
small  fraction  of  the  CIO, 

CIO  +  hv(266  nm)  -  CI  +  O  (29) 
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The  decay  of  oxygen  atoms  in  the  presence  of  excess  CIO  was  followed  by 

monitoring  the  time  dependence  of  fluorescence  signal  which  was  continuously 

excited  by  a  microwave  discharge  resonance  lamp.  The  lamp  was  operated  with 

a  low  pressure  of  helium  (<1  Torr)  containing  a  small  fraction  of  O2. 

The  ozone  storage  bulb  contained  a  mixture  of  1%  to  2%  O3  in  nitrogen, 
while  the  CI2  was  stored  neat.  These  species  were  leaked  through  needle  valves 
into  the  main  gas  flow.  Ozone  in  the  gas  flow  was  measured  in  a  35.0  cm 
absorption  cell  that  was  placed  within  a  multipass  optical  arrangement.  Using 
modified  White  cell  optics,  30  passes  of  the  254  nm  Hg  line  from  a  pen-ray  lamp 
were  sent  through  the  absorption  cell  for  an  effective  path  length  of  10.50  m. 
The  chlorine  was  measured  in  a  216  cm  absorption  cell  using  a  single  pass  of  the 
366  nm  Hg  line,  also  from  a  pen-ray  lamp.  Both  atomic  lines  were  isolated  using 
suitable  band  pass  filters.  Typically,  the  absorption  cells  were  upstream  from  the 
reaction  cells,  although  in  a  few  experiments  the  CI2  and  O3  were  measured  after 
the  flow  exited  the  reaction  cell.  Because  CI2  absorbance  at  254  nm  was  not 
totally  negligible  [a  =  1.6  x  10"^^  cm^  (Ref.  60)],  the  reference  light  intensity  for 
the  [O3]  determination  was  always  measured  with  CI2  flowing.  The  CI2 
absorption  cross  section  at  366  nm  and  the  O3  absorption  cross  section  at  254 
nm  were  taken  to  be  1.01  x  lO'^^  cm^  (Ref.  60)  and  1.147  x  lO'^'^  cm^  (Ref.  61- 
63),  respectively. 

The  Pyrex  reaction  cell  measured  16  cm  along  its  longer  axis  and  had  an 
internal  diameter  of  4  cm.  The  two  laser  beams  counterpropagated  along  the 
longer  axis.  Around  the  middle  of  the  cell  were  four  1.5  cm  diameter  side  arms. 
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each  perpendicular  to  the  long  axis  of  the  cell  and  at  90°  to  each  other.    The 

resonance  lamp  radiation  entered  the  cell  through  one  of  the  side  arms  and  the 

fluorescence  signal  was  collected  through  a  neighboring  arm.     The  central 

portion  of  the  cell  was  surrounded  by  a  jacket  through  which  thermostated 

liquids  were  flowed  to  control  the  temperature  of  the  gas  mixture  inside  the 

reactor.    The  gas  mixture  entered  the  cell  through  several  ports  very  near  the 

window  at  one  end  of  the  long  axis  and  exited  the  cell  through  simileir  ports  near 

the  opposite  window.  Because  the  chemistry  initiated  by  the  excimer  laser  beam 

completely  titrated  one  component  (O3)  of  the  gas  mixture  within  much  of  the 

cell  volume,  the  cell  was  designed  to  have  minimum  total  volume  and  low  dead 

space,  i.e.,  gas  flowed  through  all  volume  elements  of  the  cell  at  approximately 

equal  rates.   The  typical  linear  flow  rate  through  the  cell  was  14  cm  s'^  and  the 

repetition  rate  of  the  two  laser  sequence  was  usually  0.4  Hz.  Therefore,  the  gas 

mixture  within  the  entire  volume  of  the  reaction  cell  was  replenished  between 

excimer  laser  pulses.    The  temperature  of  the  gas  mixture  was  measured  by 

replacing  one  of  the  end  windows  with  an  acrylic  flange  through  which  a 

copper-constantan  thermocouple  could  be  inserted.   The  errors  in  the  reported 

temperatures  are  estimated  to  be  no  more  than  ±  1.0  K  at  the  extreme 

temperatures  and  less  at  intermediate  temperatures. 

Oxygen  resonance  lamp  radiation  was  focused  into  the  reaction  zone  by  a 

2-inch  focal  length  MgF2  lens.    The  reaction  zone  was  viewed  by  a  solar  blind 

photomultiplier  tube  through  a  similar  lens.  The  volumes  between  the  resonance 

lamp  and  the  reaction  cell,  and  between  the  reaction  cell  and  the  photomultiplier 
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tube  were  purged  with  a  mixture  of  1%  O2  in  nitrogen.  This  excluded  room  air 

and  also  acted  as  a  filter  of  extraneous  emissions  from  the  resonance  leimp.   A 

CaF2  window  between  the  cell  and  the  photomultiplier  tube  eliminated  the 

possibility  of  hydrogen  atom  detection.  Fluorescence  signals  were  accumulated 

using  photon  counting  techniques  in  conjunction  with  multichannel  scaling. 

Each  sweep  of  the  analyzer  was  triggered  simultaneously  with  the  excimer  laser. 

From  50  to  500  flashes  were  averaged  to  obtain  sufficient  signal  to  noise  ratio 

for  quantitative  kinetic  analjrsis. 

The  total  pressure  in  the  flow  system  was  measured  with  a  capacitance 
manometer.  Due  to  the  necessarily  fast  flow  rate  and  the  small  (4.0  mm  i.d.) 
tubing  connecting  the  veu-ious  components  of  the  flow  system,  there  were 
measurable  pressure  gradients  between  the  absorption  cells  and  the  reaction  cell. 
Quantitative  adjustments  were  made  for  these  gradients  in  the  calculation  of  the 
CI2  and  O3  concentrations  in  the  reaction  cell  under  each  set  of  conditions.  The 
magnitude  of  the  adjustment  was  largest  at  the  lowest  pressure  (15%  at  25  Torr) 
and  negligible  at  the  highest  (1%  at  200  Torr).  The  nitrogen  buffer  gas 
comprised  at  least  94%  of  the  mixture  for  all  experiments  and  its  flow  rate  was 
monitored  using  a  calibrated  electronic  mass  flowmeter. 

The  concentration  of  CIO,  the  excess  species  in  this  technique,  was  derived 
from  the  concentration  of  ozone  as  measured  in  situ.  Therefore,  it  was  not 
necessary  to  have  an  absolute  calibration  of  the  photolysis  laser  fluences. 
However,  it  was  important  to  know  that  following  the  excimer  laser  pulse  the 
condition  [CI]  >  [O3]  held  throughout  the  reaction  zone.  Therefore,  the  excimer 
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laser  fluence  was  measured  in  each  experiment.   As  will  be  discussed  below,  it 

was  also  important  to  monitor  the  266  nm  laser  fluence.  Both  of  these  quantities 

were  determined  using  a  photodiode-based  calibrated  radiometer.  With  the  front 

optic  of  the  excimer  laser  approximately  2  m  from  the  center  of  the  reaction  cell, 

the  beam  was  rectangular  in  cross  section  and  measured  2.0  cm  by  4.0  cm. 

When  measured  through  a  0.5  cm  diameter  aperture,  the  fluence  peaked  at  the 

center  of  the  beam  and  dropped  off  by  10%  per  0.5  cm  distance  from  the  center 

in  both  the  vertical  and  horizontal  directions.  The  beam  from  the  Nd:YAG  laser 

was  aligned  at  the  center  of  the  volume  irradiated  by  the  excimer  laser.  It  was 

estimated  to  be  0.4  ±  0.1  cm  in  diameter. 

The  reagent  purities  £ind  sources  were  as  follows:  N2  (99.999%,  Spectra 
Gases,  Inc.);  CI2  (99.9%,  Matheson  Gas  Products,  Inc.);  O2  (99.99%,  Spectra 
Gases,  Inc.).  Ozone  was  prepeired  in  a  commercial  ozonator  using  UHP  oxygen. 
It  was  stored  at  195  K  on  silica  gel  and  degassed  at  77  K  before  use.  The  other 
gases  were  used  without  further  purification. 

Results 

In  the  absence  of  competing  reactions  that  either  deplete  or  enhance  the 
ground  state  oxygen  atom  (0(^P))  concentration,  the  temporal  behavior  of 
[0(^P)]  following  the  266  nm  laser  pulse  can  be  described  by  the  relationship 

In  {[0(3p)]t/[0(3p)]t-}  =  -  (k2[C10]  +  k'd)(t  - 1')  =  -  k'(t-t')     (XVIII) 
where 

k'd  =  kaoCCy  +  ki8  (XIX) 
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In  Equation  XIX,  k3o  and  k^s  are  the  rate  coefficients  for  the  following  processes: 

0(3p)  +  CI2  -  products-  (30) 


0(3p)  -♦  loss  by  diffusion  from  the  viewing 
zone  and  reaction  with  background 
impurities  (18) 


A  typical  experimental  0(^P)  temporal  profile  is  shown  in  Figure  18. 
Unexpectedly,  a  significant  build-up  and  decay  of  0(^P)  occured  before  the  266 
nm  laser  fired;  possible  sources  of  this  0(^P)  and  its  implications  for  our  study 
of  reaction  (2)  are  discussed  below.  It  should  be  noted  that  the  vertical  axis  in 
Figure  18  has  units  of  concentration.  To  construct  Figure  18,  the  fluorescence 
signal  before  the  266  nm  laser  fired  was  scaled  to  account  for  the  fact  that  the 
351  nm  laser  photolyzed  the  entire  field  of  view  of  the  detection  system,  while 
the  266  nm  laser  photolyzed  only  15%  of  the  detector  viewing  zone.  The  size  of 
the  viewing  zone  was  estimated  by  placing  a  series  of  apertures  in  front  of  the 
351  nm  beam  and  noting  the  variation  of  fluorescence  signal  strength  with  beam 
size. 

Typical  decays  of  0(^P)  generated  by  the  266  nm  laser  pulse  are  shown  in 
Figure  19.  At  each  temperature  and  pressure  a  minimum  of  five  and  an  average 
of  eight  experiments  were  performed  at  various  values  of  [Os]^,.  Because  the 
Nd:  YAG  laser  fluence  was  measured  in  each  experiment,  the  amount  of  CIO  lost 
via  266  nm  photolysis  could  be  quantified.   The  expression 
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Figure  18.  An  experimental  0  atom  temporal  profile  obtained 
under  the  following  conditions:  [CI2]  =  1.40  x  10^^ 
molecules  cm'^;  [O^]^  =  7.27  x  10^^  molecules  cm"^;  [Cl]o 
=  1.82  X  10^'^  molecules  cm"^;  total  pressure  =  25  Torr; 
temperature  =  298  K;  multichannel  scalar  dwell  time  = 
25  microseconds. 
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Figure  19.  Typical  O  atom  temporal  profiles.  These  experiments 
were  carried  out  under  the  following  conditions  (all 
concentrations  expressed  in  molecules  cm'^):  25  Torr 
total  pressure;  T  =  298  K;  [Clg]  =  9.8  x  lO^S;  [Qgl^  = 
1.73  X  1013  (A),  5.15  X  1013  (g)^  8.76  x  10^3  (C);  and  [Cll^ 
=  1.10  x  10l4  (A),  1.66  X  IQl^  (B),  2.29  x  lO^^  (C). 
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[ClOlt  =[03]„  -  [0(3p)]t-  (XX) 

was  used  to  calculate  the  CIO  concentration  under  the  assumption  that  no  loss 
of  CIO  occurred  during  the  delay  between  photolysis  laser  pulses.  The  ratio 
[0(^P)]/[C10]  immediately  following  the  266  nm  laser  pulse  was  determined  for 
each  experiment,  and  was  tjrpicaily  in  the  range  0.01-  0.04.  Hence,  it  was 
appropriate  to  correct  each  measured  decay  rate  for  slight  deviations  from 
pseudo-first  order  conditions.  These  corrections  were  derived  from  computer 
simulations  of  reactions  (2),  (30)  and  (18)  under  the  range  of  experimental 
conditions  employed.  (A  discussion  of  the  computer  simulations  appears  in 
Appendix  A.)  Best  fit  values  for  each  decay  rate  (k')  were  obtained  from  Unear 
regression  analyses  of  the  experimental  data  over  at  least  two  1/e  times.  Each 
value  of  k'  was  then  increased  by  2%  or  less  using  the  appropriate  nonpsuedo- 
first  order  correction.  Representative  plots  of  [CIO]  versus  k'  are  shown  in 
Figure  20.  These  data  were  subjected  to  linear  least  squares  analyses  to  give 
values  for  k2.  The  results  are  presented  in  Table  VI.  Note  the  separate  columns 
for  k2 (uncorrected)  and  k2(corrected).  The  former  represents  the  best  fit  to  the 
data  when  k'  was  not  corrected  for  nonpseudo-first  order  behavior  and  [CIO]  was 
set  equal  to  [03]o  less  the  amount  photolyzed  to  produce  [0(^P)].  The  latter  k2 
values  include  the  small  nonpseudo-first  order  correction  to  k'  and  additional 
corrections  to  [CIO]  discussed  below. 

The   absorption   cross   section   vised   to   calculate   the   amount   of  CIO 
photolyzed  by  the  Nd:YAG  laser  was  estimated  experimentally.  The  signal  level 
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Figure  20.  A  k'  versus  [CIO]  plot  of  typical  data  taken  at  25  Torr 
total  pressure  of  N2  and  at  three  temperatures.  Note 
that  the  298  and  231  K  data  have  been  displaced  upward 
by  1000  and  2000  s"^,  respectively.  Solid  lines  are 
obtained  from  linear  least  squares  analyses  and  give  the 
following  rate  coefficients  (in  units  of  10"^^  cm^ 
molecule-1  s'l):  4.99  at  231  K,  4.07  at  298  K,  3.44  at  367 
K. 
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Table  VI.   Summary  of  k2  Determinations 


Temperature 
(K) 

Pressure 
(Torr) 

kg  X  lOll 
(uncorrected)  ^^^) 
(cm^molecule'^s"^) 

k2  X  lOll 

(corrected)  ^*^ 

(cm^molecule'^s'^) 

231 

25 

4.64  +  0.16 

4.99  ±  0.17 

238 

25 

4.25  ±  0.43 

4.47  ±  0.43 

252 

200 

3.10  ±  0.22 

3.97  ±  0.30 

255 

25 

4.15  ±  0.20 

4.42  ±  0.20 

255 

200 

3.25  ±  0.20 

3.95  ±  0.32 

257 

25 

3.86  ±  0.15 

4.10  +  0.11 

257 

25 

4.30  ±  0.15 

4.47  +  0.19 

275 

25 

3.94  +  0.17 

4.10  +  0.16 

298 

25 

3.55  ±  0.21 

3.84  +  0.14 

1     298 

25 

3.47  ±  0.25 

3.59  +  0.25 

298 

25 

3.51  ±  0.14 

3.62  +  0.13 

298 

25 

3.87  ±  0.36 

3.93  +  0.36 

298 

25 

3.89  ±  0.15 

4.07  +  0.17 

298 

50 

3.91  ±  0.23 

4.07  +  0.24 

298 

50 

3.62  ±  0.15 

3.98  +  0.18 

298 

50 

3.70  ±  0.12 

3.91  +  0.11<'=> 

298 

50 

3.53  ±  0.19 

3.76  +  0.14 

298 

50 

3.76  ±  0.16 

3.98  +  0.14 

298 

100 

3.73  ±  0.68 

3.97  +  0.67 

298 

200 

3.39  ±  0.27 

3.73  +  0.24 

(a)  Errors  are  2a. 

(b)  Uncorrected  values  have  not  been  adjusted  for  non  pseudo-first  order 
conditions  and  [CIO]  loss  by  reaction  (31).   See  text  for  details. 

(c)  Carried  out  with  283  nm  photolysis  of  CIO,  under  "reversed"  and  "normal" 
flow  conditions. 
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Table  VI.   (Continued) 


Temperatiire 
(K) 

Pressure 
(Torr) 

kg  X  lOll 
(uncorrected)  ^^^) 
(cm^olecule'^s"^) 

k2  X  lOll            1 
(corrected)  (*) 
(cm^molecule"^  s'  ^ ) 

298 

200 

3.39  ±  0.28 

3.71  ±  0.24 

298 

200 

3.55  +  0.12 

3.76  ±  0.13 

298 

500 

3.18  ±  0.16 

4.03  ±  0.20 

338 

25 

3.16  +  0.60 

3.26  ±  0.12 

359 

200 

2.99  ±  0.14 

3.09  +  0.15 

360 

25 

2.89  ±  0.13 

2.96  +  0.12 

1            367 

25 

3.35  ±  0.13 

3.44  ±  0.13 

(a)  Errors  are  2a. 

(b)  Uncorrected  values  have  not  been  adjusted  for  non  pseudo-first  order 
conditions  and  [CIO]  loss  by  reaction  (31).   See  text  for  deteiils. 
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immediately  after  the  Nd:YAG  laser  fired  was  directly  proportional  to  the 

concentration  of  oxygen  atoms.  If  the  Nd:  YAG  laser  was  not  preceded  by  a  pulse 
from  the  excimer  laser,  then  the  photolyte  was  O3.  If  all  the  O3  had  been 
converted  to  CIO  via  reaction  with  CI  atoms  created  in  the  351  nm  excimer 
pulse,  then  the  signal  was  due  to  CIO  photolysis.  Therefore,  in  back-to-back 
experiments  with  constant  [CI2],  [O3]  and  266  nm  fluence,  the  ratio  of  0(^P) 
signal  with  and  without  351  nm  photolysis  should  be  identical  to  the  ratio  of  the 
CIO  and  O3  absorption  cross  sections  at  the  Nd:YAG  laser  wavelength  (the 
0(^D)  product  of  O3  photolysis  is  rapidly  quenched  to  0(^P)  by  N2).  This 
experiment  resulted  in  a  value  of  0.35  for  the  signal  ratio.  If  (7(03,266  nm)  is 
taken  to  be  9.0  x  lO'^^  cm^  (Ref.  60)  then  a(C10,266  nm)  is  approximately  3.1  x 
10'^^  cm^.  This  is  somewhat  lower  than  the  low  resolution  cross  section  at  this 
wavelength  that  is  reported  in  the  literature.^  However,  we  carried  out  the 
identical  signal  level  comparison  near  283  nm,  a  wavelength  where  the  high 
resolution  cross  section  has  been  characterized.^^  Using  a  frequency  doubled, 
Nd:YAG  pumped  tunable  dye  laser  the  rotational  structure  in  the  CIO  spectrum 
was  reproduced  by  observing  the  resonance  fluorescence  signal  as  the  wavelength 
was  scanned.  At  282.95  nm,  the  peaks  of  the  R(19.5)  and  P(16.5)  lines  of  the 
^^^3/2  '  ^^^3/2  ^'^  band,  we  observed  a  factor  of  1.65  +_  0.20  more  fluorescence 
signal  when  the  excimer  laser  photolyzed  CI2  than  when  the  excimer  laser  beam 
was  blocked.  Based  on  the  dye  laser  linewidth  employed  and  the  literature 
vedues  for  the  ozone  and  CIO  cross  sections,  we  expected  a  ratio  of  1.80  ±_  0.40. 
This  result  confirms  our  value  for  a2QQ(C\0). 
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As  mentioned  above,  the  delay  between  the  two  laser  pulses,  t',  was 

adjusted  to  be  long  enough  for  reaction  (12)  to  go  to  completion.   The  value  of 

the  CI  +  O3  rate  constant  is  reasonably  well  known^^  so  an  appropriate  delay 

time  could  be  calculated.  As  a  check  the  delay  time  was  varied  until  a  constant 

signal  and  decay  rate  were  observed,  indicating  that  all  the  O3  had  been 

converted  to  CIO.    The  majority  of  experiments  were  carried  out  with  delay 

periods  of  either  3.4  ms  or  6.4  ms. 

The  temporal  behavior  of  CIO  could  be  monitored  by  following  both  the 
0(^P)  signal  level  produced  by  266  nm  photolysis  and  the  measured  value  of  k'. 
For  some  conditions  (lower  temperature,  higher  pressures)  it  was  observed  that 
CIO  was  decreasing  on  the  time  scale  of  the  delay  between  the  two  lasers.  The 
disappearance  of  CIO  is  probably  due  to  self-reaction,  a  process  that  has  not  been 
completely  characterized.^^  In  order  to  make  a  correction  for  the  amount  of  CIO 
lost  during  the  delay  between  laser  pulses  and  also  during  the  O  atom  decay,  a 
series  of  experiments  were  carried  out  in  which  the  delay  time  was  varied  at 
fixed  [O3],  [CI2]  and  laser  fluences. 

For  the  process 

CIO  +  CIO  -*  products  (31) 

the  time  dependence  of  [CIO]  is  described  by  the  relationship 

2k3it  =  [ClOlt-l  -  [ClO]^-!  (XXI) 
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If  we  define  ks^-  to  be  the  second  order  rate  coefficient  for  CIO  removal  under 

our  experimental  conditions,  then  plotting  [010]^"^  versus  time  should  jdeld  a 
straight  line  of  slope  2k3i'  and  intercept  equivalent  to  [C10]o"^  The  absolute 
concentration  of  CIO  needed  to  determine  k3i  was  deduced  in  two  manners.  In 
the  first  method,  Equation  XVII  was  rearranged  to  calculate  [CIO].  Here  k2  was 
taken  from  fixed  delay  experiments  at  the  same  experimental  temperature  and 
pressure.  k3i'  was  determined  from  Equation  XXI.  Then  the  vedues  for  [CIO] 
in  the  fixed  delay  experiments  were  adjusted  to  account  for  the  loss  of  CIO 
during  the  delay.  This  procedure  was  iterated  until  constant  values  of  k2  and 
k3i'  were  obtained.  The  second  method  involved  relating  the  signal  level 
immediately  following  the  266  nm  laser  pulse  to  the  value  of  [CIO].  The  signal 
calibration  could  be  arrived  at  by  two  procedures.  When  enough  data  was 
available  under  the  same  experimental  conditions,  a  [CIO]  versus  signal 
calibration  curve  was  constructed.  Otherwise,  signal  versus  delay  time  was 
plotted  and  extrapolated  to  time  zero  and  that  signal  level  normalized  to  [03]^. 
This  normalization  factor  was  then  used  to  put  the  signal  level  at  each  delay 
time  on  an  absolute  concentration  scale.  Because  of  the  reciprocal  relation  in 
Equation  XXI,  the  value  of  k^y  is  quite  sensitive  to  the  scaling  factor  vised  in 
setting  the  value  of  [CIO].  Also  the  signal  level  is  quite  dependent  on  the 
operating  conditions  of  the  resonance  lamp.  Therefore,  the  values  of  k3i' 
determined  from  the  first  method  are  more  precise.  However,  when  using  k'  as 
a  measure  of  [CIO],  we  assume  that  O  atoms  are  removed  only  by  reaction  with 
CIO  and  CI2;  this  assumption  appears  to  be  justified.   A  typical  plot  of  [ClO]'^ 
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versus  time  appears  in  Figure  21.     The  end  results  of  a  number  of  such 

experiments  appear  in  Table  VII.  These  results  are  presented  only  as  a  measure 

of  the  phenomenological  loss  rate  of  CIO  in  our  system,  not  as  a  definitive 

measurement  of  k3i. 

Hayman  et  al.^^  report  298  K  values  for  ksi  which  are  substantially  faster 

than  the  values  for  k3i'(298  K,P)  we  have  determined.    However,  our  255  K 

results  appear  to  agree  reasonably  well  with  the  low  temperature  values  of  ksi 

of  Hayman  et  al.  According  to  Hayman  et  al.,®^  reaction  (31)  has  an  important 

branch  to  form  a  weakly  bound  dimer  which  can  either  decompose  or  react 

rapidly  with  chlorine  atoms.   The  relevant  reaction  scheme  is  given  below: 

CIO  +  CIO  +  M  *♦  (C10)2  +  M  (31a,-31a) 

CIO  +  CIO  -*  OCIO  +  CI  (31b) 

CIO  +  CIO  ^  other  products  (31c) 

CI  +  (C10)2  -*  CI2  +  ClOO  (32) 

ClOO  +  M  -  CI  +  O2  +  M  (33) 

CI  +  ClOO  -  CI2  +  O2  (major)  (34a) 

CI  +  ClOO  -*  2C10  (minor)  (34b) 

CI  +  OCIO  -*  2C10  (35) 
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Hayman  et  al.^^  obtained  their  kinetic  data  from  a  molecular  modulation  study. 

They  observed  much  smaller  apparent  values  for  k3i  during  the  lights-off  cycle 

(no  chlorine  atoms  present)  than  during  the  lights-on  cycle  (large  concentration 

of  chlorine  atoms  present).    Hence,  a  plausible  explanation  for  the  difference 

between  our  k3i'(298  K)  determinations  and  the  k3i(298  K)  values  reported  by 

Hayman  et  al.®^  is  that  k.3ia(298  K)  was  substantially  faster  than  k32[Cl]  under 

our  experimental  conditions,  thus  facilitating  CIO  regeneration  by  dimer 

decomposition;  this  is,  of  course,  a  favorable  situation  for  measurement  of  k2 

but  would  result  in  systematic  underestimation  of  k3i.    No  clear  variation  of 

ksi'  with  [CI]  is  evident  in  our  data,  although  our  experiments  spanned  a 

rather  narrow  range  of  chlorine  atom  concentrations  (Table  VII). 

In  order  to  make  appropriate  corrections  for  the  loss  of  CIO  during  the 

delay  between  laser  pulses,  a  reaction  system  consisting  of  reactions  (12)  and 

(31)  was  modeled  under  a  variety  of  initial  conditions  using  literature  values 

for  ki2(T)^^  and  setting  ksi  =  k^y-  (A  discussion  of  how  this  reaction  sequence 

was  modeled  can  be  found  in  Appendix  A.)    From  these  calculations  a  set  of 

correction  factors  (F)  could  be  derived: 

Ft,  =  [C10]t,/[03]o  (XXII) 

Values  of  Ff  at  t'  =  3.4  ms  are  given  in  Table  VIII  as  a  function  of 
temperature,  pressure,  [03]o,  and  [CIIo/LOsIq.  Note  that  the  maximum 
correction  is  made  for  the  highest  value  of  [CIO]  for  a  given  experiment;  also, 
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the  longer  the  delay  used  in  an  experiment  the  larger  the  correction  applied. 

For  a  few  experimental  conditions  where  k^y  was  great  enough  (i.e.,  high 
pressure,  low  temperature)  a  correction  was  made  to  the  observed  k'  for  loss 
of  CIO  during  the  decay  itself;  this  correction  never  exceeded  1.4%. 

As  a  further  check  on  the  consistency  of  our  experiments,  the  rate 
coefficient  for  reaction  (2)  was  determined  at  298  K  in  50  Torr  N2  using  283 
nm  photolysis  of  CIO  rather  than  266  nm  photolysis.  The  additional  photolysis 
wavelength  was  provided  by  the  frequency  doubled,  Nd:YAG  pumped  tunable 
dye  laser.  Also,  to  ensure  that  neither  O3  nor  CI2  were  being  lost  in  the  flow 
system,  a  measurement  of  k2(298  K)  was  carried  out  with  the  two  absorption 
cells  plumbed  downstream  from  the  reaction  cell  rather  than  in  the  "normal" 
upstream  position.  Neither  of  these  variations  in  experimental  parameters 
affected  the  observed  kinetics. 

As  mentioned  above,  a  significant  O  atom  signal  was  generated  subsequent 
to  the  excimer  laser  pulse  (Figure  18).  Leu^^  and  Vanderzanden  and  Birks^^ 
have  observed  O  atoms  from  the  reaction  of  chlorine  atoms  with  ozone.  The 
following  chemistry  was  proposed  to  explain  their  observations: 

CI  +  O3  ^  CIO  +  02(bl2:g+)  (12') 

02(bl2  +>  +  O3  -  O  +  2O2  (36) 
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Subsequently,  Choo  and  Leu^^  monitored  02(b^2g'*")  directly  by  observing  its 

near  infrared  emission  and  put  an  upper  limit  of  -0.05%  on  the  02(b^2g"'")  yield 
from  reaction  (12).  Such  a  yield  is  much  too  small  to  account  for  the  observed 
levels  of  0(^P)  produced  in  the  CI  +  O3  studies. ^^,67  Choo  and  Leu  have 
suggested  that  the  O  atoms  may  be  generated  by  reactions  of  vibrationally 
excited  CIO  formed  in  reaction  (12),  i.e., 

CI  +  O3  -  C10*(v'<18)  +  O2  (12") 

followed  by 

C10*(v'>12)  +  O3  -  O2  +  O  +  CIO  (37) 

or 

C10*(v'>3)  +  CI  ^  CI2  +  O  (38) 

In  our  system  CI  atoms  are  in  excess  over  O3  and,  therefore,  in  excess  over  the 
CIO  created  in  reaction  (12).  Also,  the  occurence  of  reaction  (38)  results  in  loss 
of  two  CIO  molecules  which  otherwise  would  have  been  present  at  t'  while  the 
occurence  of  reaction  (37)  results  in  loss  of  only  one  CIO  molecule.  Hence, 
reaction  (38)  would  potentially  have  the  larger  effect  on  the  CIO  concentration 
in  our  experiments.  To  examine  the  potential  role  of  reaction  (38)  in 
consuming  CIO,  we  simulated  the  system  chemistry  using  a  Gear  routine  (see 
Appendix  A)  to  solve  the  rate  equations  numerically.  The  following  scheme  was 
modeled: 
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CI  +  Og  -  CIO*  +  O2  (12") 

CI  +  CIO*  -  CI2  +  O  (38) 

CIO*  +  M  -*  CIO  +  M  (M=N2)  (39) 

O  +  CIO  -»  O2  +  CI  (2) 

CIO  +  CIO  -  products  (31') 

O  +  CI2  ^  CIO  +  CI  (30a) 

Although  reactions  (12")  and  (38)  may  have  other  channels,  this  scheme  was 
devised  to  have  the  greatest  impact  on  the  [CIO]  at  the  end  of  the  reaction 
period;  therefore,  only  the  branches  that  deplete  CIO  were  used.  It  should  be 
noted  that  a  large  dependence  of  the  extraneous  0(^P)  signal  on  the  total 
pressure  in  the  system  was  observed.  This  observation  would  be  consistent 
with  quenching  of  either  CIO*  or  02(b^2g"'")  by  N2.  k2  and  kgo  are  known,  and 
the  other  rate  coefficients  were  adjusted  to  reproduce  the  magnitude  and 
temporal  behavior  of  the  O  atom  signal  observed.  Typical  results  are  shown  in 
Figure  22.  It  was  found  that  the  above  reaction  scheme  led  to  no  more  than 
a  two  percent  perturbation  in  the  concentration  of  CIO.  Of  course  there  is  no 
direct  proof  that  the  assumed  reaction  scheme  is  correct.     However,  it  is  the 
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Figure  22.  Comparison  of  oxygen  atom  temporal  profile  immediately 
following  351  nm  laser  pulse  to  computer  simulation  of 
reaction  scheme  discussed  in  the  text.  Experimental 
conditions  are  as  follows:  T  =  298  K;  [Oglo  =  1.65  x  lO^^ 
molecules  cm"^;  [CI2]  =  1.42  x  10^^  molecules  cm""^;  [Cl]o  = 
1.81  X  10^"*  molecules  cm""^;  pressure  =  25  Torr.  The 
reaction  rate  coefficients  used  in  the  simulation  are:  k2  = 
3.66  X  10-11;  ki2-  =  1.18  x  lO-^;  kgo  =  2.91  x  lO'l^;  k^y  = 
3.0  X  10-14;  kg^  =  2.0  X  10-^;  kgs  =  5.8  x  IQ-l^.  All  rate 
coefficients  are  in  units  of  cm"^  molecule-1  s-1. 
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worst  case  of  the  suggested  possibilities  and  any  corrections  would  be  quite 

small.  Since  the  mechanism  for  O  atom  formation  following  the  excimer  laser 

pulse   is   not  well   understood,   the   effect   of  this   chemistry  on   the   CIO 

concentration  weis  not  teiken  into  account  in  the  determination  of  "corrected" 

values  for  kg  (Table  VI). 

The  combination  of  CI2  photolysis,  CI  +  O3  reaction,  and  CIO  photolj^is 
could  result  in  some  heating  of  the  gas  in  the  reaction  zone.  Calculations  which 
assume  worst  case  conditions,  i.e.,  P  =  25  Torr  N2,  [Cllg  =  6  x  10^'*  molecules 
per  cm^,  [O3]  =  1  x  10^'*  molecules  per  cm^,  [0]^'  =  2  x  10^^  atoms  per  cm^, 
and  all  excess  energy  dissipated  as  heat,  show  that  laser  heating  of  the  reaction 
zone  could  not  have  exceeded  2  Kelvin  in  any  experiment.  This  potential 
systematic  error  is  negligibly  small  so  it  was  not  incorporated  into  the  data 
analysis. 

As  discussed  above,  several  corrections  were  made  to  either  the  observed 
k'  values  or  in  the  calculation  of  [CIO]  from  [03]o.  For  clarity  this  set  of 
corrections  is  reiterated: 

a)  A  quantitative  correction  was  made  in  k'  for  nonpseudo-first  order 
conditions  during  the  O  atom  decay.  In  only  a  few  cases  did  this  exceed 
a  1%  adjustment. 

b)  [CIO]  was  corrected  for  the  amount  of  CIO  lost  to  photolysis  at  266  nm. 
This  correction  was  dependent  upon  knowledge  of  the  laser  fluence, 
which  was  monitored  in  every  experiment  using  a  calibrated  radiometer, 
and  on  our  estimated  value  of  a266(C10).  Because  there  was  only  a  small 
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adjustment  to  [CIO]  (on  the  average  3%)  the  final  results  were  not  very 

sensitive  to  this  correction.    For  example,  in  an  experiment  where  the 

fraction  of  CIO  photolyzed  was  above  the  average,  an  increase  in 

a266(C10)  of  a  factor  of  2  was  found  to  change  the  final  value  of  k2  by 

only  4%. 

c)  Using  the  loss  rate  of  CIO  determined  in  the  same  system,  a  correction 
was  made  to  [CIO]  for  the  CIO  that  undergoes  self-reaction  (or  other  loss 
processes)  during  the  delay  time  between  laser  firings.  The  largest 
corrections  were  made  at  higher  pressures,  lower  temperatures  and  in 
long  delay  experiments.  Because  the  simulation  of  the  production  and 
loss  of  CIO  was  sensitive  to  the  errors  in  our  measurements  of  k3i',  this 
correction  has  a  rather  large  uncertainty. 

d)  A  few  decays  were  corrected  for  CIO  loss  during  the  decay  itself; 
however,  this  correction  was  insignificant  under  most  experimental 
conditions. 

e)  It  was  concluded  that  of  the  known  possibilities  for  the  source  of  O 
atoms  prior  to  CIO  photolysis,  none  could  have  had  more  than  a  two  per 
cent  effect  on  [010]^;  no  corrections  were  made  for  this  chemistry. 

Even  though  there  were  several  corrections  made  in  order  to  reach  a  final 
value  for  k2  at  each  temperature  and  pressure,  the  magnitudes  of  the 
corrections  were  small  in  most  cases  (see  Table  VI),  the  corrections  could  be 
quantitatively  applied,  and  in  general,  the  results  are  self-consistent.     An 
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Arrhenius  plot  of  our  data  appears  in  Figure  23.  An  unweighted  least  squares 

analysis  of  all  data  yields  the  expression 

kgCT)  =  (1.68  +  0.31)  X  10-11  exp  {(241  +  53)/T)}  cm^molecule-ls-l  (XXIII) 

where  the  errors  represent  2a,  precision  only,  and  Oj^  =  AcTi^y^.  For  the  25  Ton- 
data  only,  the  expression 

k2(T)  =  (1.55  +  0.33)  X  lO'^^  exp  {(263  ±  60)/T)}  cm^molecule-ls-l    (XXIV) 

is  obtained  from  an  unweighted  least  squares  analysis.  The  difference  between 
the  Arrhenius  expression  obtained  from  the  25  Torr  data  and  that  obtained 
from  the  complete  data  set  is  primarily  due  to  the  slightly  lower  rate 
coefficients  obtained  at  high  pressure  (200  Torr)  and  low  temperature  (252-255 
K).  These  rate  coefficients  required  rather  large  corrections  for  contributions 
from  reaction  (31)  and,  therefore,  are  more  likely  to  be  in  error  than  rate 
coefficients  obtained  at  higher  temperatures  and/or  lower  pressures.  For  this 
reason,  we  believe  the  Arrhenius  expression  obtained  from  the  25  Torr  data 
only  should  be  preferred.  The  absolute  accuracy  of  k2  at  any  temperature 
within  the  range  studied  is  estimated  to  be  ±_  20%. 

Discussion 
In  all  prior  investigations  reaction  (2)  was  studied  by  flow  tube  techniques 
at  pressures  less  than  10  Torr.     The  results  of  all  studies  are  summarized  in 
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Table  IX.    In  the  earliest  study  Bemand  et  al.,'*^  measured  k2(298  K)  using 

resonance  fluorescence  to  monitor  0(^P)  in  excess  CIO.  The  CIO  was  produced 

by  the  reaction 

CI  +  OCIO  -  2  CIO  (35) 

assuming  a  stoichiometric  factor  of  2.  These  workers  reported  k2(298  K)  =  5.3 
±  0.8)  X  10'^^  cm^olecule'^s"^.  They  also  measured  the  rate  coefficient  by 
following  the  decay  of  CIO  mass  spectrometricedly  in  excess  oxygen  atoms  and 
obtained  the  result  (5.7  ±  2.3)  x  10"^^  cm^olecule'^s"^.  In  a  subsequent  study 
in  the  same  laboratory  Clyne  and  Nip^^  measured  the  temperature  dependence 
of  k2.  Again,  0(^P)  was  monitored  by  resonance  fluorescence  in  excess  CIO,  the 
latter  species  being  generated  via  reaction  (12).  The  room  temperature  rate 
coefficient  was  in  good  agreement  with  their  previous  study.  They  report  a 
significant  activation  energy  and  quote  the  Arrhenius  expression  k2(T)  =  (1.07 
±  0.30)  X  10-10  exp  {-(224  ±  76)/T}  cm^molecule'ls-l. 

Zahniser  and  Kaufman^  ^  measured  the  temperature  dependence  of  the 
ratio  k2/ki2-  Using  a  value  of  ki2  measured  directly  in  the  same  system  these 
workers  report  k2(T)  =  (3.38  ±  0.50)  x  lO"!!  exp  {(75  ±  40)/T}  cm^olecule"! 
s-1. 

The  next  reported  investigation  of  k2(T)  was  performed  by  Leu^^  using 
resonance  fluorescence  detection  of  0(^P)  in  excess  CIO.  The  CIO  radicals  were 
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produced    using    three    different    source    reactions    in    order    to    validate 

stoichiometric  assumptions  necessary  to  arrive  at  CIO  concentration  levels. 

Reactions  (12),  (35)  and 

CI  +  CI2O  -  CIO  +  CI2  (40) 

were  the  three  sources  used  at  room  temperature.  Reaction  (12)  was  used  at 
all  other  temperatures.  Leu's  value  for  k2(298  K)  is  lower  than  the  previously 
reported  values  and  he  measured  a  small  positive  activation  energy  with  k2(T) 
=  (5.0  +  1.0)  X  10-11  exp  {-(96  ±  20)/T}  cm^molecule-ls'l.  The  use  of  reaction 
(12)  as  the  only  CIO  source  in  all  experiments  at  T  ^  298  K  could  result  in  a 
systematic  error  in  Leu's  reported  temperature  dependence.  At  the  low 
pressures  employed  in  Leu's  study,  the  CIO*  +  CI  reaction  may  have  competed 
favorably  with  CIO*  deactivation,  thus  leading  to  overestimation  of  [CIO]. 
Since  the  ratio  of  k38/k39  may  be  temperature  dependent,  such  an  effect  could 
have  been  more  important  at  one  end  of  the  investigated  temperature  range 
than  at  the  other  end.  It  should  be  emphasized  that  while  the  abovementioned 
systematic  error  in  Leu's  [CIO]  determination  is  possible,  there  currently  exists 
insufficient  information  concerning  CIO*  chemistry  to  prove  or  disprove  this 
conjecture. 

Schwab  et  al.^^  employed  an  experimental  apparatus  in  which  both 
reactants  (CIO  by  laser  magnetic  resonance  and  0(^P)  by  resonance 
fluorescence)  and  one  product  (CI  atoms  by  resonance  fluorescence)  could  be 


110 
monitored.    [Interestingly,  this  is  the  only  O  +  CIO  study  where  CIO  was 

directly  measured  in  the  reaction  zone.]   Again,  a  somewhat  lower  value  was 

measured  (k2(298  K)  =  (3.5  +  0.5)  x  10'^^  cm^olecule-^s-^)  and  these  workers 

observed  essentially  no  temperature  dependence  for  k2. 

Ongstad  and  Birks^'^^  measured  k2(T)  in  a  discharge  flow  system  using 

the  same  three  sources  of  CIO   as  Leu^^.      0(^P)   was   followed  via  the 

chemiluminescence  from  NO2*  generated  by  reacting  the  oxygen  atoms  with 

NO  added  to  the  detection  region  of  their  flow  tube.  These  workers  measured 

k2(T)  directly  and  also  relative  to  the  reaction 

O  +  NO2  -  NO  +  O2  (41) 

In  a  successive  measurement  scheme  k4i(T)  was  also  measured.  The  relative 
measurements  yielded  somewhat  higher  values  presumably  due  to  nonpure 
source  gases  (CI2O,  CIO2,  O3)  or  other  channels  for  the  source  reactions. 
Ongstad  and  Birks  measured  a  value  for  k2(298  K)  that  agrees  with  the  other 
more  recent  studies  and  a  small  "negative  activation  energy."  They  reported 
the  expression  k2(T)  =  (2.61  ±  0.60)  x  lO'^l  exp  {(97  +  64)/T}  cm^olecule"! 
s-1. 

The  only  other  study  of  k2  at  or  near  room  temperature  reported  in  the 
literature  was  by  Margitan.^^  CIO  radicals  were  generated  in  a  flow  tube  via 
reaction  (40).  Downstream  from  this  source  0(^P)  was  created  by  laser 
photolysis  of  the  CIO  and  followed  by  resonance  fluorescence.     The  CIO 
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concentration  was  measured  directly  in  the  flow  by  absorption.  However,  large 
corrections  (up  to  20%)  had  to  be  made  for  CIO  loss  between  the  CIO  detection 
region  and  the  0(^P)  detection  region.  Margitan  used  literature  values  for  k3i 
to  make  these  corrections.  Given  the  recent  advances  in  our  understanding  of 
reaction  (31),^^  a  large  uncertainty  mxist  be  associated  with  the  magnitude  of 
Margitan' s  correction  for  CIO  loss  via  the  self-reaction.  His  results  are  also 
very  dependent  on  the  value  chosen  for  the  CIO  absorption  cross  section. 
Margitan  reports  that  E/R  lies  within  the  range  +  200  K  and  k2(298  K)  =  (4.2 
+  0.8)  X  10-11  cm^olecule-ls-1. 

As  seen  by  the  comparison  in  Table  IX,  there  is  very  little  difference 
between  the  value  of  k2(298  K)  from  our  experiments  and  from  any  of  the  other 
recent  studies.  However,  at  the  lower  temperatures  typical  of  the  middle 
stratosphere  our  results  indicate  significantly  faster  values  for  k2(T)  than  any 
of  the  other  recent  investigations  (see  Table  IX).  Hence,  model  calculations 
which  employ  our  expression  for  k2(T)  would  predict  somewhat  larger  ozone 
depletion  due  to  chlorofluorocarbon  injection  than  calculations  which  take  k2(T) 
from  previously  available  data. 

Summary 

We  have  measured  k2  as  a  function  of  temperature  and  pressure.  Our 
results  indicate  a  lack  of  any  pressure  dependence  at  298  K  over  the  range  25 
Torr  to  500  Torr.  Although  our  298  K  rate  coefficient  agrees  well  with 
previous  studies,  our  observation  of  an  activation  energy  that  is  more  negative 
than  any  previously  reported  leads  to  a  significant  difference  between  our  result 
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and  other  recent  measurements  at  temperatures  relevant  to  stratospheric 

chemistry. 


CHAPTER  V 
CONCLUSIONS 

In  the  preceding  chapters,  our  studies  of  the  temperature  dependence 
of  the  rate  coefficients  for  the  reactions 

O  +  HO2  -  OH  +  O2  (1) 

and 

O  +  CIO  -  CI  +  O2  (2) 

are  presented.  Each  chapter  contains  a  comparison  of  our  measurements 
with  all  previously  published  work  on  these  reactions.  Subsequent  to 
publication  of  our  results,  there  have  been  no  additional  measurements  of 
these  rate  coefficients.  However,  it  is  useful  to  discuss  the  current 
"evaluated"  rate  coefficient  for  each  reaction,  as  well  as  recent 
developments  relevant  to  our  studies. 

As  discussed  in  Chapter  III,  there  is  very  good  agreement  on  the 
room  temperature  value  of  k^  among  the  recent  direct  measure- 
ments. ^^'^'''^^'^^'^^  Also,  the  temperature  dependence  for  k^  observed  in 
our  study  is  essentially  identical  with  that  seen  by  Keyser.^®  Therefore, 
the  current  NASA  evaluation^^  of  ki(T)  is  an  average  of  the  five  studies 
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at  room  temperature  fitted  to  the  temperature  dependence  from  our  work 

and  that  of  Keyser.^^ 

In  order  to  extract  a  value  for  k^CT)  from  our  data,  small  corrections 

were  made  for  some  secondary  chemistry.  The  set  of  reactions  that  were 

included  in  our  modeling  of  the  secondary  processes  is  given  in  Table  III 

in  Chapter  III.   The  values  for  the  rate  coefficients  listed  in  Table  III  are 

very  nearly  identical  with  the  current  NASA  evaluations^^  except  for  the 

rate  coefficient  for  reaction  (21). 

OH  +  HO2  -*  H2O  +  O2  (21) 

A  recent  study  by  Keyser^^  demonstrated  that  previous  low  pressure 
studies  of  reaction  (21)  in  discharge  flow  experiments  had  been  subject  to 
secondary  chemistry  problems  leading  to  low  values  for  the  rate 
coefficients  (see  Ref.  69  for  discussion).  The  current  NASA^^ 
recommended  value  at  room  temperature  is  approximately  50%  faster 
than  the  value  of  k^  used  in  our  study.  A  re-evaluation  of  the  effect  of  the 
secondary  chemistry  listed  in  Table  III  but  with  the  current  value  for  k2i 
leads  to  no  more  than  a  few  percent  change  in  the  worse  case  for  any  rate 
coefficient  determination;  therefore,  it  is  not  necessary  to  make  any 
changes  in  our  reported  results.  The  methods  used  to  simulate  the 
secondary  chemistry  in  these  studies  is  discussed  in  the  final  section  of 
Appendix  A. 
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For  reaction  (2),  the  NASA  evaluation^®  is  based  on  a  least  squares 

analysis  of  the  seven  most  recent  measurements  of  k2.^^'^^"^^    There  is 

reasonable  agreement  for  k2(298  K)  but  some  uncertainty  still  exists  as  to 

the  magnitude  and  direction  of  the  temperature  dependence. 

In  the  0(^P)  +  CIO  study  small  corrections  were  made  for  the 
amount  of  CIO  lost  to  266  nm  photolysis.  These  corrections  were  based 
on  an  absorption  cross  section  value  that  is  somewhat  lower  than  that 
given  in  a  recent  evaluation^®  of  the  literature  on  the  CIO  absorption 
spectrum.  However,  as  was  pointed  out  in  Chapter  IV,  a  factor  of  2 
increase  in  the  value  of  (7(C10,266  nm)  changes  k2  by  only  4  percent. 
Recently,  the  CIO  cross  section  at  266  nm,  near  the  peak  of  the 
continuum,  has  been  found  to  be  independent  of  temperature.^® 
Temperature  independence  had  been  assumed  in  our  study. 

During  our  study  of  reaction  (2)  (see  Chapter  IV),  we  were  able  to 
estimate  values  for  the  process 

CIO  +  CIO  -  product  (31) 

under  different  conditions  of  temperature  and  pressure.  These  results  are 
presented  in  Table  VII.  When  our  study  was  completed,  only  the  work  of 
Hayman  et  al.^^  was  available  for  comparison.  Subsequently,  two  more 
measurements  of  ksi  have  appeared  in  the  literature. ^^'^^  The  results  of 
each  of  these  studies  of  k3i(P,T)  are  compared  with  our  reported  values 
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in  Table  X.  The  reasonably  good  agreement  between  our  reported  values 

for  the  phenomenological  loss  of  CIO  and  the  recent  literature  values  for 
k3i  give  support  to  the  interpretation  of  the  CIO  loss  in  our  experiments. 
In  summary  it  can  be  stated  that  two  of  the  main  goals  of  this 
research  have  been  reached.  The  temperature  dependence  of  the  rate 
coefficients  for  the  atmospherically  important  reactions  (1)  and  (2)  have 
been  measured.  Our  values  for  k^CT)  and  k2(T)  are  currently  included  in 
the  primary  chemical  kinetic  database  for  use  in  stratospheric  modeling. 
Secondly,  we  have  clearly  demonstrated  that  laser  flash  photolysis- 
resonance  fluorescence  is  a  viable  technique  for  use  in  kinetics  studies  of 
radical-radical  reactions.  Also,  it  has  been  demonstrated  that  neither  of 
these  radical-radical  reactions  are  pressure  dependent.  Even  though  the 
entrance  channel  for  the  reactions  are  attractive,  the  possible  bound 
species,  HO3*  and  CIO2*,  are  not  stabilized  on  the  time  scale  of  collisions 
at  pressures  approaching  one  atmosphere. 
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APPENDIX  A 
FURTHER  DETAILS  OF  THE  EXPERIMENTAL  METHODS 

Introduction 
One  of  the  shortcomings  of  most  journal  articles  in  experimental 
physical  chemistry  is  that  space  only  permits  a  cursory  treatment  of  the 
experimental  method  employed  in  the  work.  A  reader  interested  in  more 
details  is  left  with  options  such  as  searching  out  any  references  to  similar 
experiments,  generalized  descriptions  of  the  method  in  textbooks  or 
review  articles,  or  establishing  direct  contact  with  the  researchers. 
Because  the  work  described  in  the  preceding  three  chapters  represents  a 
somewhat  innovative  approach  to  studies  of  radical-radical  reactions, 
more  information  was  given  concerning  the  experimental  method 
(especially  in  Chapter  II)  than  would  be  found  in  a  published  account  of 
a  more  standard  flash  photolysis  resonance  fluorescence  experiment. 
However,  many  details  were  excluded.  This  appendix  contains  additional 
information  about  some  aspects  of  the  experiments  and  should  redress 
some  of  the  deficiencies  of  the  experimental  method  sections  of  the 
preceding  three  chapters. 
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Flow  Systems 

The  components  of  the  various  gaseous  reaction  mixtures  were  leaked 

into  the  flow  system  through  needle  valves  and  monitored  using  calibrated 

electronic  mass  flowmeters  (Matheson  Gas  Products,  Inc.).  Using  1/4"  o.d. 

teflon  tubing  with  either  Swagelock  or   Cajon   fittings,   the   separate 

components  were  combined  and  flowed  through  a  sequence  of  cells  (the 

reaction  cell  and  UV  absorption  cells).    The  gas  mixture  was  exhausted 

through  a  rotary  vane  mechanical  pump  (Sargent  Welch  Model  1402  or 

1397).  The  pressure  of  the  system  and  flow  rate  through  the  reaction  ceil 

was  controlled  by  adjusting  the  flow  of  each  component  into  the  system 

and  by  use  of  a  throttle  valve  on  the  exhaust  port.    A  separate,  larger 

exhaust  valve  could  be  used  to  bypass  the  throttle  valve  to  facilitate  leak 

testing  and  flushing  of  the  system.  A  Pyrex  trap  held  at  the  boiling  point 

of  nitrogen  (77  K)  was  plumbed  directly  above  the  vacuum  pump.    This 

prevented  corrosive  condensibles  (H2O2,  O3,  CI2)  from  reaching  the  pump 

and  degrading  the  pump  oil.      The  cold  trap  also   prevented  volatile 

components  of  the  pump  oil  from  backstreaming  into  the  system  under 

pump  down.     It  is  worth  noting  that  ozone  has  a  significant  vapor 

pressure  at  77K  and  therefore  cannot  be  completely  trapped.     O3  is 

particularly  damaging  to  red  rubber  vacuum  hose,  with  small  perforations 

appearing  in  the  tubing  even  after  limited  exposure.  However,  the  O3  can 

be  destroyed  in  a  length  of  heated  copper  tubing  plumbed  upstream  of  any 

hose  and  the  pump. 
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All  experiments  were  carried  out  under  "slow  flow"  conditions.  This 

prevented  the  depletion  of  any  components  of  the  reaction  mixture  and 

the  build-up  of  reaction  products  that  might  become  involved  in  secondary 

chemical    processes.      The    experiments    were    conducted    in    a   signal 

averaging  mode,  i.e.,  signal  traces  following  individual,  sequential  laser 

pulses  were  accumulated  to  improve  the  counting  statistics.    Therefore, 

the  lower  limit  of  the  flow  rate  was  set  by  the  condition  that  every  laser 

flash  encountered  a  fresh  gas  mixture.   Because  the  photolysis  was  down 

the  long  axis,  the  entire  volume  of  the  cell  had  to  be  replenished  between 

laser  shots.  This  was  especially  critical  in  the  case  of  the  O  +  CIO  study. 

As  explained  in  Chapter  IV,  one  component  of  the  gas  mixture,  O3,  was 

completely  titrated  after  each  laser  pulse.    Therefore,  for  the  O  +  CIO 

experiment,  a  low  volume,  low  dead  space  cell  was  designed  that  would 

ensure  total  flushing  of  the  cell  volume  in  a  convenient  time  period.  The 

linear  flow  rate  was  typically  14  cm/s  and  the  laser  repetition  rate  was  0.4 

cm/s  in  the  O  +  CIO  experiments.   In  the  O  +  HO2  studies,  only  a  small 

fraction  of  the  gas  mixture  was  photolyzed  so  that  the  condition  of 

complete  gas  exchange  between  laser  pulses  was  not  as  stringent. 

The  flowmeters  were  calibrated  during  each  study  using  one  of  two 

different  techniques.  In  the  first  technique,  the  pressure  drop  in  a  known 

volume  was  timed  while  monitoring  the  flowmeter  reading  generated  by 

a  fixed  flow  out  of  that  volume.   In  the  second  method,  the  gas  output  of 

the  flowmeter  was  passed  through  a  buret  fitted  with  a  latex  reservoir 
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filled  with  soap  solution.  Introduction  of  a  small  amount  of  the  soap  into 

the  gas  flow  created  thin  soap  films  across  the  buret  that  could  be  timed 

as  the  'bubble'  rose  through  the  buret.    It  is  worth  pointing  out  that  as 

long  as  all  the  flowmeters  used  in  an  experiment  were  calibrated  near  the 

same   time   with  the   same  method,   then   the  meters   should  be   well 

calibrated  relative  to  each  other.  Calibration  on  an  absolute  scale  was  not 

necessary  because  it  was  the  fractional  flow  through  each  flowmeter  that 

was  used  in  any  concentration  determination.    In  none  of  these  studies 

were  the  measured  rate  coefficients  directly  dependent  on  the  calibration 

of  the  flowmeters.    In  the  O  +  HO2  studies  the  [H2O2]  was  measured  in 

situ  by  UV  photometry.    The  flowmeter  readings  were  used  to  calculate 

total  flow  rate  and  the  value  of  [O3]  in  each  experiment.  In  the  O  +  CIO 

study  both  CI2  and  O3  were  monitored  via  UV  photometry.    In  fact,  for 

most  experiments  neither  the  CI2  flow  nor  the  flow  of  the  O3/N2  mixture 

were  even  measured  by  flowmeters.     These  flows  were  only  a  small 

fraction  of  the  total  flow  and  more  constant  concentrations  of  CI2  and  O3 

were  possible  when  the  concentrated  O3  flow  and  pure  CI2  flow  were  not 

passed  through  the  heated  and  restrictive  stainless  steel  flowmeters. 

Direct  measurement  of  the  critical  low  concentration  components  in 

the  flow  greatly  reduces  the  uncertainty  in  any  measured  rate  coefficient. 

Any    systematic    loss    of   these    components    due    to    dark    reactions, 

heterogeneous  processes,  or  thermal  decomposition  will  be  detected  and 

taken  into  account.  Of  course,  accurate  absorption  coefficients  are  then 
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important.  In  the  O  +  HO2  studies,  however,  it  is  actually  the  ratio  of  the 

H2O2  absorption  coefficients  at  the  monitoring  wavelength  and  at  the 
photolysis  wavelength  that  is  critical.  If  the  wavelengths  are  reasonably 
close  and  the  cross  sections  are  taken  from  a  single  study,  then  the  ratio 
of  cross  sections  can  be  more  accurate  than  either  individual  cross  section. 
Because  a  titration  was  used  in  the  O  +  CIO  study,  the  rate  coefficient 
was  only  directly  dependent  upon  the  concentration  of  a  single  species 
and,  therefore,  a  single  absorption  cross  section,  i.e.,  that  of  O3  at  253.7 
nm. 

The  pressure  in  the  system  was  measured  with  an  absolute 
capacitance  manometer  (MKS  Instruments  Model  220).  Because  small 
tubing  connected  the  larger  components  of  the  flow  system  and  the  flow 
rate  was  relatively  large,  a  pressure  gradient  existed  between  the 
absorption  cell(s)  and  the  reaction  cell.  Therefore,  the  concentration  of 
the  species  measured  by  photometry  were  corrected  by  the  appropriate 
ratio  of  pressures.  The  pressure  ratios  were  measured  by  alternately 
coupling  the  pressure  gauge  directly  to  the  reaction  cell  and  the 
absorption  cell  for  each  set  of  pressure  and  flow  conditions  used  in  the 
experiments. 

Resonance  Lamps 

The  lamp  itself  is  essentially  the  same  as  described  by  Davis  and 
Braun.*^^  A  1"  diameter  MgF2  window  is  attached  to  one  end  of  the  lamp 
using  Torr  Seal  low  vapor  pressure  epoxy  (Varian  Associates).  Compared 
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to  MgF2,  LiF  is  more  transmitting  in  the  vacuum  UV,  but  it  is  a  much 

weaker  material.  Helium  containing  a  small  amount  of  O2,  either  as  an 
impurity  or  from  a  separate  reservoir,  is  flowed  rapidly  through  the  lamp 
at  low  pressure,  a  few  Torr.  Helium  is  used  rather  than  argon  in  order 
to  avoid  the  extraneous  emissions  of  the  latter.  For  these  experiments, 
a  rotary  vane  pump  (Sargent  Welch  Model  1397)  pulled  the  helium 
through  the  lamp  and  a  LN2  cooled  trap.  A  100  watt,  2.45  GHz 
microwave  power  supply  (KIVA  Instruments  Model  MPG4)  was  used  to 
maintain  the  discharge  in  an  Evenson  cavity^'*  (Ophthos  Instruments). 
For  atomic  species,  the  gas  flow  through  the  discharge  should  always  be 
in  the  direction  away  from  the  window  in  order  to  lessen  reabsorption  by 
ground  state  atoms.  Emissions  by  other  species  can  be  reduced  by  using 
higher  purity  gases  or  possibly  filtering  the  lamp  emissions.  The  quality 
of  the  resonance  fluorescence  signal  is  very  dependent  upon  the  operation 
of  the  lamp.  The  concentration  of  O2  in  the  helium,  the  total  pressure, 
the  power  level  of  the  discharge,  the  position  of  the  discharge  along  the 
axis  of  the  lamp  and  the  tuning  of  the  cavity  are  all  factors  in  determining 
the  detection  sensitivity  of  O  atoms.  The  optimum  set  of  conditions  can 
only  be  found  by  systematically  varying  all  these  parameters  and 
empirically  arriving  at  the  best  signal  level.  While  there  is  a  temperature 
difference  maintained  between  the  emitting  atoms  in  the  lamp  and  the 
absorbing  atoms  in  the  reaction  cell,  the  fluorescence  will  still  be 
proportional  to  the  concentration  of  atoms  in  the  cell  as  long  as  the 
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emission  profile  of  the  lamp  does  not  change  during  the  measurement.  By 

merely  replacing  oxygen  with  another  appropriate  precursor,  the  lamps 

can  be  used  as  sources  of  resonance  radiation  for  any  number  of  elements. 

Two   final  points   are   noted   on  the  operation  of  these   lamps.      The 

transmission  of  the  window  may  degrade  under  use  due  to  production  of 

trapped  absorbers  in  the  crystal.    Also,  care  must  be  taken  to  keep  all 

exposed  surfaces  clean,  since  almost  any  adsorbed  contaminate  will  have 

a  vacuum  UV  absorption. 

Fluorescence  Detection 

The  fluorescence  from  the  O  atoms  in  the  reaction  cell  excited  by  the 

resonance  lamp  was  detected  using  a  solar  blind  photomultiplier  tube 

(either  a  Thorn  EMI  Model  G26L314LF  or  EMR  Schlumberger  Model 

5181ER  or  542G).  The  photocathode  material  in  these  end-on  tubes  is  Csl 

coated  on  the  inside  of  the  MgF2  window.    The  detector  was  mounted 

perpendicular  to  both  the  path  of  the  photolysis  laser  and  the  direction 

of  light  from  the  resonance  lamp.  Because  the  resonance  wavelengths  are 

in  the  vacuum  UV  (130-135  nm),  MgF2  or  CaF2  optics  were  employed  to 

couple  light  into  and  out  of  the  reaction  cell.    The  space  between  the 

resonance  lamp  output  window  and  the  reaction  cell  as  well  as  the  space 

between  the  cell  and  the  window  of  the  photomultiplier  tube  were  purged 

with  a  gas  transparent  to  the  resonance  wavelengths.  The  purge  gas,  N2, 

excluded  room  air,  which  contains  O2,  H2O  and  other  absorbers.   Helium 

should  never  be  used  as  a  component  of  the  purge  gas  in  front  of  a 
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photomultiplier  tube,  because  it  can  permeate  the  window,  or  window 

seals,  hence,  destroying  the  tube.    A  small  controlled  amount  of  a  filter 

gas,  O2,  could  be  added  to  one  or  both  of  the  purges  to  increase  the  ratio 

of  the  resonance  radiation  to  the  background  emissions  from  the  lamp. 

Because  light  was  coupled  into  and  out  of  the  reaction  cell  using  MgF2 

lenses  (ca.  2"  focal  length),  the  distances  of  the  resonance  lamp  discharge 

and  the  photoactive  area  of  the  PM  tube  from  the  reaction  zone  were  also 

critical  parameters  in  the  detection  of  O  atoms. 

Signal  Acquisition 

Data  were  accumulated  in  a  fashion  usually  referred  to  as  photon 
counting.  Weak  pulses  from  the  photomultiplier  tube  were  fed  into  a 
photon  discriminator  (Mechtronics  Nuclear  Model  551).  To  bias  against 
the  weaker  signal  arising  mostly  from  thermal  processes,  only  pulses 
larger  than  an  adjustable  threshold  level  triggered  the  discriminator  to 
output  a  TTL  logic  pulse.  These  logic  pulses  were  the  input  for  the  signal 
analyzer  (either  a  Tracor  Northern  TN-1500  or  NS-570)  operating  in  the 
multichannel  scaling  (MCS)  mode.  The  analyzer  accumulates  the  data  in 
a  time  versus  signal  counts  grid.  The  time  axis  is  subdivided  into  equal 
increments  sized  to  match  the  rate  of  the  decay  process  observed. 

An  external  trigger  signal  marked  the  start  of  each  MCS  sweep. 
Multiple  sweeps  were  accumulated  in  order  to  improve  the  signal 
statistics.  The  scheme  used  to  synchronize  the  MCS  sweep  with  the  laser 
pulses  was  different  in  each  of  the  studies.     In  the  first  O   +   HO2 
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experiment,  a  Lumonics  Model  TE-261-2  was  used  as  the  photolysis  laser. 

Timing  pulses  from  a  pulse  generator  in  the  laser  power  supply  were 
routed  to  a  digital  delay  generator  (California  Avionics  Laboratories 
Model  lOlAR).  The  delay  generator  put  out  two  pulses,  a  t^  pulse  that 
was  practically  synchronous  with  the  faster  pulse  from  the  laser  triggered 
the  signal  analyzer  and  a  delayed  pulse  triggered  the  firing  of  the  laser. 
The  delay  between  these  pulses  was  typically  10  ms.  Triggering  the 
analyzer  prior  to  the  laser  flash  allows  determination  of  the  background 
light  prior  to  generation  of  O  atom  signal.  In  the  second  O  +  HO2 
experiment,  a  Lambda-Physik  EMG200  excimer  laser  was  employed.  A 
separate  pulse  generator  (Global  Specialties  Model  4001)  was  used  for  the 
master  pulse.  This  triggered  the  digital  delay  generator  which  again  put 
out  a  Iq  pulse  to  the  signal  analyzer  and  a  delayed  pulse  to  the  laser.  In 
the  O  +  CIO  experiment,  two  lasers  were  used  but  a  similar  scheme  was 
employed.  A  pulse  generator  integral  to  the  power  supply  of  the  Nd:YAG 
laser  (Quantel  481A)  supplied  the  master  pulse  to  trigger  the  delay 
generator.  In  this  case,  the  delay  generator  t^  pulse  triggered  both  the 
Lambda-Physik  excimer  laser  and  the  signal  analyzer.  A  delayed  pulse 
then  triggered  the  flash  lamps  of  the  YAG  laser. 

After  sufficient  signal  had  been  accumulated,  the  background  was 
subtracted  in  the  signal  analyzer  and  the  data  downloaded.  In  the  first 
two  studies,  a  printout  of  signal  counts  versus  channel  number  was 
obtained.    As  a  first  cut  analysis,  this  data  was  hand  plotted  on  semi- 
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logarithmic  graph  paper.   Because  the  signal  decays  were  expected  to  be 

exponential  then  a  plot  of  the  natural  logrithm  of  the  signal  versus  time 

yields  a  line  whose  slope  is  equal  to  the  decay  fate  for  that  experiment. 

Later,  more  precise  analysis  was  performed  by  entering  the  data  into  a 

personal  computer  for  linear  least  squares  analysis.     In  the  O  +  CIO 

experiment,  an  interface  was  devised  that  allowed  direct  transfer  of  the 

data   from   the   TN-1500   signal   analyzer   into   a   Radio   Shack   Color 

Computer  II  for  subsequent  analysis. 

UV  Photometry 

In  all  the  experiments  described  in  the  preceding  three  chapters,  the 

concentrations  of  the  minor,  but  critical,  constituents  in  the  gas  flow  were 

monitored  by  UV  photometry.    As  discussed  above,  the  accuracy  of  the 

measured  rate  coefficient  is  determined  in  part  by  the  accuracy  in  the 

concentration  of  critical  components  in  the  gas  mixture.  In  the  O  +  HO2 

studies,  [H2O2]  was  monitored  directly  in  the  gas  flow  via  its  absorption 

at  either  202.6  nm  or  228.8  nm.    The  concentration  of  ozone,  the  other 

minor  component  of  the  gas  mixture  in  the  O  +  HO2  experiments,  was 

calculated  based  on  the  total  pressure  in  the  system  and  the  fractional 

flow  from  a  dilute  O3/N2  mixture.  However,  the  fractional  concentration 

of  O3  in  the  12  liter  Pyrex  storage  bulb  containing  the  ozone  mixture  was 

periodically  determined  by  measuring  the  absorption  at  253.7  nm  due  to 

a  known  pressure  of  the  mixture.    For  the  O  +  CIO  study,  both  CI2  and 
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O3  were   monitored   in  situ  via  absorption  at  366  nm  and  254  nm, 

respectively. 

The  absorption  cells  were  Pyrex  tubing  with  size  25  0-ring  joints  as 
end  pieces  and  1/4"  o.d.  side  tubes  near  each  end  to  allow  for  the  gas  flow, 
UV  quality  quartz  windows  were  clamped  across  each  of  the  0-ring  joints. 
Depending  on  the  wavelength  to  be  monitored,  three  different  lamps  were 
used.  A  zinc  hollow  cathode  lamp  was  the  source  of  a  zinc  ion  line  at 
202.6  nm.  Because  the  output  of  the  zinc  lamp  includes  several  other 
strong  lines  in  the  same  wavelength  region  it  was  necessary  to  isolate  the 
202.6  nm  line  with  a  monochromator.  Typically,  the  lamp  was  mounted 
6"  from  a  1  cm^  aperture  on  one  window  of  the  absorption  cell.  The 
entrance  slit  of  a  1/4  m  Jarrell-Ash  monochromator  was  positioned  near 
the  other  window  of  the  cell.  With  the  slits  used  in  these  measurements, 
the  monochromator  was  able  to  completely  resolve  the  202.6  nm  Zn"*"  line 
from  the  nearest  line  at  206.2  nm  (also  Zn"*").  A  quartz  envelope 
photomultiplier  tube  (RCA  1P128  or  similar)  in  a  light  tight  box  was 
mounted  directly  to  the  monochromator  at  the  exit  slit.  The  high  voltage 
(typically  700  V  d.c.)  was  divided  equally  down  the  dynode  chain  and  the 
current  measured  with  an  electrometer.  For  a  larger  absorption  cell,  100 
cm  or  more,  a  few  tenths  of  a  microamp  of  signal  could  be  expected  for  the 
unattenuated  light  level. 

At  the  other  photometric  wavelengths,  the  light  sources  were 
brighter  and  the  spectra  less  congested  near  the  analytical  wavelengths. 
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Therefore,  the  desired  lines  could  be  isolated  using  commercially  available 

bandpass  filters.  A  mercury  pen-ray  lamp  (UV  Photonics)  was  used  to 
monitor  ozone  at  254  nm  and  also  for  chlorine  at  366  nm.  The  former  line 
dominates  the  lamp  output  and  care  should  be  taken  to  avoid  having  too 
much  light  entering  the  absorption  cell.  Distance  and  black  vinyl  masks 
are  useful  means  of  attenuation.  The  bandpass  filters  were  obtained  from 
Acton  Research  Corporation.  Both  filters  had  stated  maximum 
bandwidths  of  12  nm  (full  width  at  half  maximum).  The  254  nm  filter  is 
approximately  10%  transmitting  at  the  peak,  while  the  365  nm  filter  is 
25%  transmitting.  The  bandpass  filters  were  mounted  directly  over  the 
opening  in  the  box  containing  the  PM  tube  and  this  assembly  positioned 
near  the  exit  window  of  the  absorption  cell.  Typically,  more  signal  at 
lower  voltages  could  be  expected  for  detection  of  the  Hg  atomic  lines  than 
for  the  zinc  line  described  above.  The  lamp  output  at  366  nm  is  actually 
three  separate  lines  at  365.0  nm,  365.5  nm,  and  366.3  nm.  The  CI2 
spectrum  is  quite  broad  in  this  region,  therefore,  resolution  of  this  triplet 
feature  is  not  necessary. 

In  the  O  -I-  CIO  study,  it  was  necessary  to  have  quite  a  long  path 
length  for  the  measurement  of  ozone.  This  was  accomplished  by  passing 
the  254  nm  output  of  the  Hg  lamp  30  times  through  a  35  cm  long 
absorption  cell.  The  multipassing  was  done  using  modified  White  cell 
optics. ^^  The  White  cell  mirrors  had  a  broad  band  reflective  coating  that 
included  254  nm  and  the  windows  of  the  cell  anti-reflection  coated  at  254 
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nm.  In  order  to  accommodate  the  large  number  of  passes,  a  wide  aperture 

was  needed.    Therefore,  the  cell  was  made  with  size  50  0-ring  joints  on 

each  end. 

There  is  a  common  interference  associated  with  the  use  of  the  254 
nm  line  from  the  mercury  pen-ray  lamp.  The  output  also  contains  a 
strong  emission  at  185  nm,  which  is  attenuated  by  O2  in  room  air. 
Photolysis  of  O2  leads  to  a  production  of  O3  which  then  attenuates  the 
254  nm  line.  The  interference  can  be  eliminated  by  constantly 
replenishing  the  air  between  the  lamp  and  the  absorption  cell,  i.e., 
directing  a  small  fan  toward  this  region. 

For  the  actinometry  experiments  carried  out  in  conjunction  with  the 
O  +  HO2  studies,  the  ozone  concentration  was  determined  in  much  the 
same  fashion  as  described  above.  However,  because  the  emission  was 
monitored  on  short  times  scales,  it  was  more  convenient  to  operate  the 
lamp  on  a  direct  current  (dc)  power  supply  rather  than  the  usual 
alternating  current  (ac)  transformer.  The  lamp  was  started  using  the  arc 
from  a  Tesla  coil  instead  of  the  high  voltage  pulse  provided  by  the 
standard  ac  supply.  Under  dc  operation,  the  lamp  intensity  eventually 
diminished  due  to  distillation  of  the  Hg  to  one  side  of  the  lamp. 
Periodically,  the  measurement  was  halted,  the  polarity  reversed,  and  the 
lamp  restarted. 
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Numerical  Simulations 

Numerical  integrations  were  used  at  several  junctures  in  these 
studies  to  assess  and  quantify  interferences  in  the  rate  coefficient 
measurement  due  to  secondary  chemistry.  It  is  worthwhile  to  give  a  brief 
description  of  these  simulation  techniques.  Because  the  chemistry  to  be 
simulated  was  relatively  simple  and  more  sophisticated  programs  were  not 
available  in  our  laboratory  at  the  time  of  the  initial  studies,  a  simple 
brute  force  approach  was  taken  in  writing  the  integration  programs. 
Initially,  a  differential  equation  was  written  to  express  the  change  in 
concentration  with  time  for  each  species  (dXj/dt).  This  set  of  equations 
was  integrated  by  repeatedly  calculating  the  change  in  each  concentration 
(AXj)  over  small  time  steps  (At),  correcting  each  concentration  for  the 
change  (Xj  +  AX{),  and  calculating  the  concentration  change  for  the  next 
time  step.  As  long  as  the  time  step  is  small  compared  to  the  rate  at  which 
the  chemistry  is  occurring,  the  integration  will  be  accurate. 

More  elegant  programming  contains  an  algorithm  for  adjusting  the 
time  steps  to  match  the  rates  as  the  integration  proceeds,  making  for  a 
more  efficient  program.  However,  for  the  relatively  simple  reaction 
schemes  discussed  in  the  preceding  chapters,  the  running  time  of  the 
integration  on  a  personal  computer  using  the  simple  approach  was  not  a 
major  constraint.  In  each  case,  the  appropriateness  of  the  time  step  was 
checked  by  demonstrating  that  a  smaller  time  step  gave  the  same  result. 
A  discussion  of  a  re-evaluation  of  the  secondary  chemistry  effecting  the 
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measurement  of  the  O  +  HO2  rate  coefficient  appears  in  the  final  chapter. 

In  this  case,  the  integration  was  performed  with  the  program  QCMP022 
purchased  from  the  Quantum  Chemistry  Program  Exchange  at  Indiana 
University.  The  agreement  between  the  results  of  this  latter  integration 
with  the  results  of  simpler  calculations  discussed  in  Chapters  II  and  III 
confirms  the  validity  of  our  initial  conclusions. 


APPENDIX  B 
LIST  OF  REACTIONS 

For  the  convenience  of  the  reader  a  list  of  the  chemical  reactions  in 
numerical  order  is  given  in  Table  XI. 
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Table  XI.   List  of  Reactions 


1.  0(3p)  +  HO2  -  OH  +  O2 

2.  0(3p)  +  CIO  -»  CI  +  O2 

3.  O2  +  hv  -  2  O 

M 

4.  O  +  O2  -*  O3 

5.  O3  +  hv  -  O  +  O2 

6.  O  +  O3  -•  2  O2 

7.  OH  +  O3  ^  HO2  +  O2 

8.  HO2  +  O3  -  OH  +  2  O2 

9.  CI  +  CH4  -  HCl  +  CH3 

10.  OH  +  HO2  -*  H2O  +  O2 

11.  H2O  +  0(lD)  -  2  OH 

12.  CI  +  O3  -»  CIO  +  O2 

13.  H2O2  +  hv  -*  2  OH 

14.  O3  +  hv  -p->  0(lD)+02(alAg) 

-»  0(3p)+02(X32^ 

15.  OH  +  H2O2  -  HO2  +  H2O 

16.  0(lD)  +  N2  ^  0(3?)  +  N2 

17.  0(3p)  +  H202-T-^  OH  +  HO2 

->  H2O  +  O2 

18.  0(3p)  -  bckgrd  loss 

19.  OgCa^Ag)  +  O3  -  0(3?)  +  2  O2 

20.  0(lD)  +  H2O2  -  OH  +  HO2 

M 

21.  OH  +  HO2  -*  H2O  +  O2 


22.  HO2  +  HO2  ■*  H2O2  +  O2 

23.  0(3?)  +  OH  -  O2  +  H 

24.  H  +  HO2  -  2  OH 

M 

25.  H  +  O2  -♦  HO2 

26.  H  +  O3  -  OH  +  O2 

27.  HO2  -•  bckgrd  loss 

28.  CI2  +  hv  -  2  CI 

29.  CIO  +  hv  -*  CI  +  O 

30.  0(3?)  +  CI2  -»  products 

31.  CIO  +  CIO  -  products 

32.  CI  +  (C10)2  -  CI2  +  ClOO 

M 

33.  ClOO  -  CI  +  O2 

34.  CI  +  ClOO-T-s.  CI2  +  O2 

->2C10 

35.  CI  +  OCIO  -  2  CIO 

36.  02(biZ+)  +  O3  ^  O  +  2  O2 

37.  C10*(v'>12)  +  O3  - 

O2  +  O  +  CIO 

38.  C10*(v'>3)  +  CI  ^  CI2  +  O 


M 

39. 

CIO*  -  CIO 

40. 

CI  +  CI2O  -  CIO  +  CI2 

41. 

0  +  NO2  -  NO  +  O2 
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